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ABSTRACT 


Synoptic  and  cl inatoloqical  dynamic  studies  generally  rely 
on  bulk  aerodynamic  flux  foraulae  to  describe  air  sea  heat  and 
■omentum  exchanqe  on  synoptic  and  cl iaatological  scales. 
Baroaetric  pressure  maps  (which  involve  an  intrinsic  temporal 
averaging  of  the  wind)  and  wind  roses  provide  two  sources  of 
spatial  and  temporal  wind  inforaation  for  flux  calculations. 
Several  investigators  have  shown  that,  due  to  the  non-linear 
dependence  of  the  bulk  aerodynamic  formulae  on  the  winds,  time- 
averaged  estimates  of  the  fluxes  based  on  vector  averaged  winds 
systematically  underestimate  the  actual  time-averaged  fluxes. 

Using  10  tc  21  years  of  three-hourly  sampled  sea  surface 
aetecrol  observations  from  9  weatherstations  in  the  North 
Atlantic  ocean  and  2  weatherstations  in  the  North  Pacific  Ocean, 
the  three-hourly  stresses,  latent  heat  fluxes  and  sensible  heat 
fluxes  were  calculated.  The  sampled  data  and  the  calculated 
fluxes  were  then  averaqed  over  periods  varying  up  to  28  days. 
The  estimates  cf  the  averaged  fluxes  based  on  the  vector 
averaqed  winds  were  then  compared  to  the  directly  averaged 
values.  -*  •" 

A  simple  analysis  revealed  that  an  upper  bound  for  the 
difference  in  the  two  stress  calculations  was  directly 
proportional  to  the  sum  of  the  x  and  y  component  wind  variances 
lost  throuqh  the  averaging  process  (in  agreement  with  Fofonoff, 
1960)  and  inversely  proportional  *o  the  square  of  the  vector 
averaqed  wind  speed.  The  wind  averaged  and  directly  averaged 
flux  estimates  were  grouped  according  to  the  Beaufort  wind  speed 
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category  and  the  period  over  which  the  variates  were  averaged.  A 
■altivariate  regression  was  then  performed  to  optimise  a 
transformation  from  the  wind  averaged  to  the  directly  averaged 
case. 

For  all  fluxes,  the  transfor aation  draaatically  improved 
the  wind  averaged  estiaates  of  the  cliaatological  leans  and 
variances  of  the  directly  averaged  fluxes.  The  residual  error 
between  the  two  estiaates  was  decreased  up  to  a  factor  of  5  over 
the  uncorrected  case  and  the  correlation  coefficients  showed  a 
moderate  increase.  The  regression  coefficients  showed  siailar 
values  for  all  teaperate  latitude  stations. 

Cased  on  consistencies  observed  in  the  wind  speed  and 
averaging  period  dependencies  of  che  multivariate  coefficients, 
an  empirical  formula  was  found  which  interpolated  the  wind  speed 
and  averaging  dependence  and  duplicated  the  aultivariate 
regression  results.  The  data  from  the  ten  temperate  latitude 
staticns  were  grouped  and  a  single  formula  found  which  cnly 
moderately  increased  the  errors  between  the  wind-averaged  and 
directly  averaged  estimates.  The  geographically  averaged  formula 
was  not  applicable  at  Station  N ,  located  at  the  northern 
extremity  of  the  Forth  Pacific  Trade  Hind  region. 

Analysis  of  the  28  day  wind-averaged  flux  spectral 
estimates  showed  that  they  underestimated  the  28  day  directly 
averaged  flux  spectral  estiaates.  Application  of  the  specific 
ship  empirical  formula  greatly  improved  agreement  between  the 
two  spectral  densities  and  reduced  the  residual  series  power 
density  at  all  frequencies.  High  latent  heat  flux  errors  at 
Station  N,  could  be  reduced  by  application  of  a  seasonal 
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The  data  were  also  grouped  into  monthly  wind  rose 
configurations  and  the  Mind  rcse  monthly  flux  estimates  Mere 
compared  to  the  directly  calculated  long-term  monthly  mean 
fluxes.  In  all  cases,  the  Hind  rose  fluxes  compared  favourably 
nith  the  directly  calculated  fluxes. 
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The  spatial  and  teaporal  scales  of  the  turbulent 
environment  for  air-sea  coupled  systeas  nay  vary  over 
siqnificant  orders  of  aaqnitude.  For  exaaple,  treads  m  the 
long-ten  heatinq  and  coolinq  of  the  Earth's  surface  of 
consequence  to  nan  aay  require  many  decades  to  deter nine  while 
turbulent  diffusion  of  aoaentua  within  the  sea  Bust  be  measured 
in  fractions  of  a  second. 

Fiedler  and  Panofsky  (  1970)  have  suggested  several 
overlapping  teaporal  (and  spatial)  scale  ranges  to  aid  in  the 
description  of  tiae-dependent  neteoroloqica 1  phenomena.  The 
first  is  the  microscale  ranging  from  less  than  a  second  to  about 
cne  hour.  The  mesoscale  ranges  from  several  minutes  to  several 
days,  the  synoptic  scales  start  at  about  one  day  and  extend  to 
several  weeks.  The  seasonal  and  climatic  scales  are  longer  than 
the  synoptic  scales. 

Hind  stress  and  sensible  and  latent  heat  fluxes  can  be 
described  as  the  vertical  transfer  of  horizontal  momentum, 
temperature,  and  humidity  from  the  air  to  tne  sea.  Microscale 
experimental  techniques  are  required  to  measure  them.  Aaong  the 
techniques  developed  are  the  eddy-correlation,  profile  and 
dissipation  methods.  A  detailed  discussion  of  them  may  be  found 
in  the  turbulence  literature.  Briefly,  however,  the  vertical 
and  downstream  wind  components  and  temperature  and  humidity 
fluctuations  are  sampled  at  a  high  rate  (usually  several  Hertz) 
for  periods  up  to  one  hour.  The  transport  fluxes  are  determined 
through  correlation  of  the  vertical  velocity  fluctuations  with 


the  downstream  wind  component,  temperature  and  humidity 
rluctuations  as  shown  below: 


T  =  - 


He  =  E 


p  v/v'  Vi' 

\ 

1.  1  (a) 

O  Cp  wt' 

1  1 

1.  1  (b) 

E  w'Q' 

1.  1(c) 

where  1  is  tne  microscale  wind  stress,  is  the  air  density 
,  hs  is  the  sensible  heat  flux.  Hi  is  the  latent  heat  flux,  w' 
is  tne  vertical  velocity,  u*  is  the  downstream  velocity 
fluctuation,  t*  is  the  absolute  air  temperature  fluctuation,  g' 
is  the  absolute  hunidity  fluctuation,  S  is  the  latent  heat  of 
vaporization  (2.46x10*  J  kg  -*)  ,  Cp  is  the  specific  heat  of  air 
at  constant  pressure  <1.0x10J  J  kg-*  ®C  ** )  and  the  overbar 
represents  the  ensemnle  average  where  all  units  are  assumed  to 
oe  s.l.  Normally  ergodicity  is  invoked  and  in  practice  a  time 
average  is  actually  used. 

Eguaticns  1.1  are  exact.  The  rapid  sampling  reguired, 
however,  (about  10  Hertz  in  the  atmosphere)  generates  too  much 
data  for  long  time  scale  studies.  Measuring  the  vertical 
fluctuation  in  particular  and  the  rapidity  of  measurements  in 
general  necessitates  using  highly  sophisticated  eguipment  that 
is  generally  jeyond  the  scope  of  the  synoptic  scale 
investigator. 

The  microscale  fluxes  can  be  related  to  the  average  of 
mesoscale  measurements  of  one  hour  (or  less)  of  sampling 


duration  through  the  bulk  parameter izations  given  by  (eg.  Roll, 


1965) 


T=  p  a  V2 


1.  2(a) 


or  in  components: 


(T*,^  "  ^ 

Ms  =  ^  Ct  V  AT 

H<_  =  Cc^L  V  A  Q 


1.2(b) 


1.2(c) 


1.2(d) 


where,  T  is  the  wind  stress  magnitude,  Tx  and  Ty  are  the  x  and  y 
components  of  the  wind  stress,  V  the  mesoscale  vector-averaged 
wind  speed,  u  and  v  are  the  x  and  y  components  of  the  mesoscale 
averaged  wind  velocity,  A  T  is  the  average  air-sea  temperature 
difference,  A  Q  is  the  average  air-sea  absolute  humidity 
difference.  Cd,  Ct  and  Cq  are  nondimensional  transfer 
coefficients  ccmmcnly  known  as  the  drag  coefficient,  and  Stanton 
and  Dalton  numbers  respectively.  They  relate  the  mesoscale  mean 
varianles  of  Equation  1.2  to  the  microscale  fluctuations  of 
Equation  1. 1. 

Except  for  specific  turbulence  studies,  hourly  wind 
measurements  are  seldom  available.  Many  calculations  of 
transport  fluxes  are  derived  from  quantities  that  are 
effectively  averaged  over  periods  much  longer  than  one  day.  If 
we  let  overbars  denote  a  tine  average  then  an  estimate  of  the 
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average  x  compete11-  stress  is: 

I 

T*.  -  x  ^  ^  *-  u  t  (.  uN  vrLtV)  1.3 

where  there  are  I  observations  and  the  subscript  i  denotes  the 
set  of  wind  components,  air  densities  and  drag  coefficient 
observations. 

An  estimate  of  the  x  component  of  stress  based  on  the 
component  averaged  winds  is  given  by: 

T*  =  C  r  £  ^i)  (_  r  £  Cck^  (.  r.  £  o.c)  * 

i(r  E4  ♦(U^YV1’' 

t 

It  is  evident  that  Eguation  1.3  and  1.4  are  not  identical.  By 
similar  arguments  the  averaged  heat  fluxes  calculated  from 
Equations  1.2(c)  and  (d)  will  differ  from  the  estimates  of  the 
averaged  heat  fluxes  based  on  the  product  of  the  averaged 
constituent  variables. 

Malkus  (  19  62)  examined  59  three-hourly  sampled  wind 
ooservations  from  the  Caribbean  Sea  and  found  that  the  fluxes 
(inexuding  momentum,  latent  heat,  and  sensible  heat)  calculated 
from  the  averaged  constituents  (Equation  1.4)  underestimated  the 
directly  averaged  fluxes  (Equation  1.3)  by  only  7. OS.  Similar 
calculations  for  the  stresses,  performed  on  three-hourly  sampled 
winds  at  Weatherstation  C  revealed  that  the  averaged  constituent 
variaoles  predicted  a  stress  magnitude  of  only  35%  of  the 
directly  calculated  value.  She  attributed  the  differences  in 
reduction  to  the  steadiness  of  the  Trade  Winds  where  the  long¬ 


term  climatolcgica 1  averages 


are  indicative  of  the  short-term 


mesoscale  and  synoptic  scale  conditions.  Over  middle  latitude 
oceans,  the  winds  are  much  more  variable,  being  dominated  by 
four  to  seven  day  cyclonic  storms,  and  the  long-term  wind 
averages  do  not  adequately  indicate  shorter  period  activity. 

The  wind  stress  input  for  many  oceanographic  calculations 
is  derived  from  estimates  of  the  geostrophic  surface  winds 
calculated  from  averaged  barometric  pressure  maps.  These  maps 
are  constructed  through  large  spatial  scale  samples  taken 
several  times  daily  then  temporally  averaged  over  periods  of  up 
to  a  month  cr  more. 

Fofonof f  (  1960)  in  calculating  oceanic  transports  in  the 
North  Pacific  Ocean  based  on  geostrophic  winds  derived  from 
monthly  pressure  maps  suggests  that  the  difference  between  the 
uiontnly  stress  and  the  stress  calculated  from  averaged  pressures 
is  proportional  to  the  variance  of  the  pressure  gradient.  The 
variance  of  the  pressure  gradient  is  equivalent  to  the  sum  of 
the  u  and  v  wind  component  variances  over  the  period  in  which 
the  pressure  was  averaged. 

Aa gaard  ( 1 S70)  provided  a  particularly  striking  example  of 
tne  effect  of  using  monthly  averaged  pressure  maps  in 
calculating  transports  in  the  Norwegian  and  Greenland  Seas. 
Starting  with  six-hourly  averaged  pressure  maps,  he  calculated 
tne  six-hourly  Sverdrup  transports  following  Fof onof f (  1960) . 
The  six-hourly  transports  were  then  averaged  over  February  1965 
to  arrive  at  a  monthly  average  transport  which  appears  in  Figure 
1(a).  He  then  averaged  the  barometric  pressures  over  the  same 
period  and  recalculated  a  pressure  averaged  estimate  of  the 
total  transport  for  February  1965.  This  result  appears  in 


Figure  1  (a) 


Figure  1  (b) 


Figure  1.  The  net  transport  in  the  Greenland  and  Norwegian  Seas 
for  February,  1965  based  on  six-hourly  barometric  pressure  maps. 
Figure  1  (a)  shows  the  month  average  of  the  six-hourly  -transports 
and  Figure  1(b)  shows  the  transports  based  on  the  month  average 
of  six-hourly  barometric  pressures.  Transport  units  are  10®  m3 
sec 'i.  (from  Aagaard,  1970) 


Figure  1(b).  It  is  guite  evident  that  in  regions  of  large 
average  transport,  the  mean  pressure  approach  underestimates  the 
mean  transport  approach  by  a  factor  of  four  to  five. 

Willebrand (  1  978)  used  twelve-hourly  averaged  pressure  maps 
to  calculate  spatial  and  temporal  wind  stress  spectra  on  the 
North  Atlantic  and  North  Pacific  Oceans.  Although  the  overall 
spectral  shapes  may  be  qualitatively  correct,  the  stress  power 
density  estimates  are  probably  low. 

Many  oceanographic  phenomena  such  as  the  reversal  of  the 
current  regime  in  the  Indian  Ocean,  the  onset  of  El  Nino  in  the 
Pdcific  Ocean,  the  reversal  of  the  California  and  Davidson 
currents  off  the  Western  United  States,  and  the  onset  of  coastal 
upwelling  have  been  linked  with  the  temporal  dependence  of  the 
wind  stress  forcing.  A  point  of  contention  is  whether  the  large 


scale  (cf  the  order  of  100  km)  wind  field  and  consequently 
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Sverdrup  transports  (curl  of  the  wind  stress)  are  important  or 
whether  local  wind  stresses  (about  30  km)  and  Ekman  transports 
are  controlling  mechanisms.  In  the  case  of  El  Nino,  it  has  been 
demonstrated  that  the  total  oceanic  stress  is  important  (see 
McCreary,  1976  and  Hurlburt,  et.al.  1976)  while  the  local  Ekman 
transports  are  important  for  coastal  upwelling  (Allen,  1980). 

For  coastal  phenomena  the  local  winds,  sampled  at  about 
hourly  intervals,  can  be  readily  obtained  from  coastal  stations. 
Open  ocean  wind  field  data  are  much  more  difficult  to  obtain 
(O'Brien,  1971)  and  geostrophic  wind  calculations  from 
oarcmetric  pressure  maps  at  present  are  the  only  viable  source 
of  wind  information  covering  long  time  and  large  spatial  scales. 
In  the  future  satellite  scatter cmet ry  may  present  an  alternative 
source  of  wind  information  but  at  present  these  measurements  are 
still  in  their  infancy. 

It  is  cnly  recently  that  time-dependent  meteorological 
conditions  have  been  used  to  drive  simulation  models  of  the  flow 
in  large  oceanic  basins.  Among  others,  a  recent  attempt  in  this 
vein  has  been  the  work  of  Huang  (1978  and  1  979)  in  modelling 
climatological  and  seasonal  variations  of  the  North  Pacific 
Ocean.  Two  stages  of  his  model  have  been  presented.  First,  the 
model  was  spun  up  tc  quasi-steady  state  using  an  initial  set  of 
data  based  on  climatolog ically  averaged  (over  20  years)  stress, 
latent  heat  flux  and  sensible  heat  flux.  The  stress  field 
resemDied  that  cf  Hellerman  ( 1967) .  In  the  next  stage,  a  second 
set  of  data  was  used  based  on  monthly  estimates  of  the  three 
fluxes  at  each  grid  point  from  monthly  averaged  values  of  air 
temperature,  vapcur  pressure  and  zonal  and  meridional  winds.  A 
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quasi-deterministic  annual  cycle  was  then  calculated  at  each 
grid  point  by  fitting  the  first  three  haraonics  of  the  annual 
cycle  to  the  second  data  set  and  accepting  the  clinatological 
information  of  the  first  data  set  for  the  constant  tern.  It  is 
evident  that  his  technique  should  lead  to  the  errors  outlined  by 
Aagaard  (  1970)  shown  in  Figure  1.  Some  of  the  inherent  averaging 
errors  in  Huang's  model  have  been  compensated  by  using  a  drag 
coefficient  of  2.5x10 -3  which  is  approximately  double  the  more 
recent  estimates  {eg.  Large  1979,  Smith  and  Banke  1974). 

Thus  an  objective  of  this  thesis  is  to  quantify  on  the 
mesoscale  and  synoptic  scale,  the  reductions  in  air  sea  fluxes 
occurring  when  estimates  are  made  using  averaged  wind  and 
temperature  data. 

Fissel  (1975)  made  an  initial  attempt  at  quantifying  the 
discrepancy  using  as  a  data  base  ten  years  of  meteorological 
observations  at  Weatherstation  P.  He  found  that  the  stress 
calculated  through  Equation  1.4  consistently  underestimated  that 
calculated  through  Equation  1.3  and  that  the  reductions  were 
quite  regular  for  averaging  periods  up  to  two  years. 

The  present  study  examines  a  data  base  of  meteorological 
observations  sampled  at  three-hourly  intervals  from  nine 
weatherstations  on  the  Atlantic  Ocean  and  two  weatherstations  on 
the  Pacific  Ocean  to  determine  the  extent  of  the  discrepancies 
referred  to  above  and  any  generalities  which  may  be  projected 
from  the  results. 

Normally  climatological  wind  stress  information  is  based  on 
a  wind  rose  analysis  as  outlined  by  Hellerman (1965) •  Here 
irregularly  sampled  wind  measurements  are  sorted  into  discrete 


wind  speed  and  direction  categories.  The  probability 


distribution  of  the  wind  coaponents  can  be  calculated  and  used 
to  obtain  estiaates  of  the  aonthly  or  annual  wind  stress. 
Calculations  of  this  sort  have  been  used  particularly  for  the 
steady  state  moaentun  flux  inputs  of  large-scale  ocean 
circulation  nuaerical  models. 

The  data  base  of  the  eleven  weatherships •  meteorological 
observations  could  be  readily  sorted  into  various  choices  of 
wind  rose  configurations.  The  “wind  rose"  fluxes  were  then 
compared  to  the  actual  individually  determined  fluxes  (using 
Equati  on  1.3)  to  investigate  the  error  induced  by  intrinsic  wind 
rose  averaging  of  the  winds  and  thus  errors  in  the  estiaates  of 
the  cliaatological  wind  stress  fields. 
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CHAPTER  II  DATA  PREPARATION 


2.1  Data  Origin 

The  data  base  for  the  study  originated  fro*  aeteorological 
observations  at  Ocean  Weatherstations  Alpha,  Bravo,  Charlie, 
Delta,  Echo,  India,  Juliet,  Kilo,  tlite,  November  and  Papa  which 
will  henceforth  be  referred  to  by  the  capital  letter  only  of  the 
station  (eg.  Weatherstation  E) .  The  quantities  recorded  at 
each  station  were  wind  speed  and  direction,  barometric  pressure, 
wet  and  dry  bulb  teaperatures  and  sea  surface  teaperatnre  which 
were  sampled  at  three-hourly  intervals.  Maps  of  the  weathership 
locations  appear  in  Figure  2. 

Data  from  Stations  A  through  H  were  transferred  froa  IBM 
punch  cards  to  seven-track  magnetic  tape  at  the  University  of 
Southaapton,  England  while  Station  N  data  were  transferred  to 
seven  track  magnetic  tape  at  Oregon  State  University.  The  wet 
and  dry  bulb  teaperatures  were  converted  to  a  dew  point 
temperature  on  these  tapes.  The  Station  P  data  were  entirely 
processed  at  the  University  of  British  Columbia  and  details  of 
the  editing  can  be  found  in  Hertzman,  et.  al.  (1974).  A  resume 
of  the  nuaber  of  years  of  data,  country  responsible  for 
collection,  and  latitude  and  longitude  of  the  stations  can  be 
found  in  Table  I.  Further  details  of  data  collection  techniques 
can  be  found  by  addressing  inquiries  to  the  national  weather 
bureau  services  cf  the  responsible  countries. 

2.2  Data  Verif icat i on 

In  order  for  tapes  froa  the  University  of  Southaapton  and 
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TABLE  I 


Details  of  the  stations  locations  and  nuaber  of  years  of  data 
that  were  collected.  The  start  date  is  1  January  of  the  start 
year  and  the  end  date  is  31  Deceaber  of  the  end  year. 


STATION)  COUNTBY  J  LATITUDE 

I  I 

)  RESPONSIBLE  | 

- + - 1 - 

A  (Netherlands(61-63  <>N 
l-U.S.  | 

B  (Canada  ( 55. 5-57.5  «N 

1-0. S.  ) 

C  |U.S.  |51.  5-53.5  <>f 

D  |U.S.  143-45  on 

£  |U.S.  134-36  on 

I  | U. K.  J59-61  on 


|U.S.  |51. 5-53.5  °N 

|0.S.  143-45  on 

|U.S.  134-36  on 

|  U. K.  J59-61  on 

|Netherlands|52.  3-54.  3  “N 

1-  U.S.  | 

I  Netherlands ) 44-46  °N 
1-Prance  | 

INorway  (65-67 

|U.S.  (28.5-31.5  oH 

(Canada  148.5-51.5  <>N 


LONGITUDE 


31-35  on 

49.  5-52.5  on 

34-37  on 
40-42  oy 
47-49  ON 
18-21  on 
17.8-20.8  on 

15-17  on 

0-4  oe 
138-142  on 
142.6-147.6  on 


STABTi  END 


DATE 

1949 

1949 

1949 

1949 

1949 

1949 

1949 

1949 

1949 

1949 

1958 


DATE 


1967 

1967 

1967 

1964 

1965 


1963 

1969 

1967 


Oregon  State  University  to  be  coapatible  with  the  University  of 
British  Colunfcia  conputer  systeas,  the  original  data  tapes  had 
to  be  reforaatted  froa  seven  to  nine  trades.  Siaultaneously, 
aissing  data  were  identified  and  flagged  by  replacing  the 
aissing  values  with  -999.9.  The  reconstructed  tine  series  were 
then  compared,  pcint-by-point,  to  extreme  values  as  outlined  in 
Table  II.  Any  data  which  had  values  outside  these  ranges  were 
considered  to  be  in  error  and  these  were  also  eliainated  from 
calculations  ny  flaqging  thea  with  -999.9.  An  example  of  the 
dew  point  temperature  time  series  for  Weatherstation  A  at  this 
level  of  processing  is  showa  in  Figure  3(a).  At  this  stage,  the 
u  and  v  wind  velocity  component  time  series  were  produced  froa 
the  wind  direction  and  speed  choosing  the  positive  x  and  y 
directions  as  east  and  north  respectively. 
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TABLE  H 

Out-of-range  Units  for  the  initial  stage  of  data  processing. 


i - 1 - 1 — - 1 - i 


I  QUANTITY 

1 

!  UNITS  J 
{  i 

LOWER 

LIBIT 

| UPPER 
t  LIBIT 

1  Hind  Direction 

“I  + 

1 0  fro»  H  | 

-180 

j  180 

1 

(Hind  Speed 

1 

|  Ba rometric 
| Pressure 

J 

(Air 

| Temperature 

1 

JDew  Point 
| Temperature 

1 

|Sea  Surface 

i  i 

1  m/sec  | 

0 

I  70 

1  1 

) millibars! 

1  1 

900 

J  (060 

1 

1  1 

1°C.  | 

1  1 

i  i 

-20 

1 

|  40 

J 

1 

1  1 

1  OC.  J 

1  1 

i  1 

-20 

)  40 

1 

I 

1  J 

l°c.  J 

-10 

|  40 

| Temperature 

1 

1  I 

1  1 

_r  L 

1 

1 

j _ 

Because  the  data  analysis  required  a  continuous  three- 
hourly  tine  series,  the  flagged  data  of  all  channels  were 
replaced  by  a  value  which  was  linearly  interpolated  between  the 
nearest  preceding  good  point  and  the  next  anteceding  good  point. 
The  wind  u  and  v  components  were  interpolated  and  the  speed  and 
direction  then  calculated  fron  the  components.  Figure  3(b) 
shows  the  interpolated  time  series  reconstructed  fron  Pigure 
3(a).  There  were  two  lain  types  of  errors  present.  The 
blackened  areas  of  Figure  3(a)  correspond  to  every  second  or 
third  measurement  being  either  missing  or  out-of -range.  Hhen 
these  values  are  replaced,  the  corresponding  sections  in  Figure 
3(b)  are  virtually  ind ist inguish  ,e  from  the  remainder  of  the 
time  series.  The  arrow  in  Figure  3(a)  indicates  where  data  are 
aissinq  for  seven  days  corresponding  to  instrument  malfunction 
or  the  ship  being  off  station.  As  can  be  seen  in  Figure  3(b), 
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the  linear  interpolation  results  in  a  flattened  region  where  the 
variance  has  been  markedly  reduced. 

The  arrows  in  Figure  3(b)  illustrate  examples  of  "spikes" 
cr  regions  where  the  data  are  discontinuous  and  obviously 
erroneous  even  though  within  the  linits  of  Table  II.  In  the 
sore  promiaant  case,  the  dew  teaperature  rises  froa  -4.8®C  to 
36.  IOC  and  back  to  -3.9<>c  in  twelve  hours.  Consequently  all  the 
data  were  visually  aonitored  on  the  U8C  Department  of 
Oceanoqraphy  PDP  -  12  coaputer.  Obvious  spikes  were  identified 
and  the  points  replaced  by  interpolated  values.  At  this  stage 
an  eight  bit  flag  recorded  each  data  point  which  had  been 
interpolated.  Table  III  qives  the  total  nuaber  of  aissing, 
erroneous,  and  out-of-range  data  points  for  each  of  the  aeasured 
quantities.  The  Station  p  inforaation  was  extracted  froa 
Fissel  (1975)  . 

With  the  exception  of  Weatherstation  I,  the  dew  point 
teaperature  displays  the  largest  number  of  interpolated  points. 
This  is  representative  of  a  large  number  of  aissing  readings 
rather  *-han  spikes.  Weathership  A ,  for  example,  has  dew  point 
measurements  for  ev«ry  second  reading  only  for  the  period  1952 
to  1954.  isiaiiar  situations  exist  for  the  other  weathersbips. 

Analyses  attempted  in  the  study  required  yearly  blocks  of 
data  that  closely  approximated  the  actual  times  series.  The 
regularly  missing  (that  is,  every  second  or  third  point)  data 
presented  little  problem  because  major  features  such  as  the 
diurnal  cycle  remained  intact.  Larqe  blocks  of  aissing  data 
severely  disserted  rhe  tinip  series  to  the  extent  that  it  was 


felt  that  they  cculd  influence  the  tinal  results. 


Por  any 
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TABLE  III 

Iotal  number  of  missing,  cut-of-range,  and  erroneous  (spiked) 
data  points  per  time  series  and  weatherstation.  The  total 
column  signifies  the  total  number  of  samples  taken  (ie.  the 
numrer  of  days  studied  times  8), 
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variate  which  contained  more  than  125  consecutive  missing  data 
points  in  any  given  month,  the  entire  year  of  data  for  that 
variaole  was  rejected.  Table  IV  shows  the  years  of  rejected 
data  for  each  ship, 

2. 3  Density  Calculations 

Page  F-S  of  the  Handbook  of  Physics  and  Chemistry  (59th 
edition)  gives  an  empirical  formula  for  the  density  of  moist  air 
as: 


P 


=  1.  23  29  (273.  13/Ia)  (P-0.  3783e/760) 


2.  1 


where  p  is  the  air  density  in  kg/m3,  Ta  is  the  dry  air 
temperature  in  °K,  P  is  the  barometric  pressure  in  mmHg,  and  e 
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The  year  (s)  of  eliminated  data  for  each  weatherstation. 


r— - 

1  SHIP 

1  A 

— *1 - 

1 

1 

-  i 

STRESS 

_r_ 

1 

1 

_  i 

HEAT 

FLUXES 

1 

1 

1 

1 

1950 

T 

1 

1 

1950 

1 

1  B 

1 

- 

1 

52,53 

1  c 

1 

- 

1 

1953 

1  o 

1 

- 

) 

52,53 

I  E 

1 

- 

i 

52,53 

|  I 

1 

1952 

1 

1952 

1  J 

1 

1952 

1 

52,57 

1  K 

1 

1953 

1 

49,50,51 

1 

i 

53,54 

1  H 

1 

1949 

1 

1949 

!  N 

1 

- 

1 

52,53 

1  P 

1 

• 

- 

l 

1 

J 

— A._ 

— a 

is  the  vapour  pressure  of  the  aoisture  in  aailg*  Bertzaan 
et.  al.  (1974)  qive  a  formula  to  calculate  the  absolute  humidity 
based  on  values  in  the  CRC  Handbook  of  Physics  and  Chemistry 
(51st  ed  ition)  : 


g=f. 4038x10*  exp  (-5  1 07.  4/Td)  2.2 

whece  q  is  the  absolute  humidity  in  gm/a3  and  Td  is  the  dew 
point  temperature  in  °K.  By  *he  ideal  gas  law; 

e/Ta=qR  2.3 

where  F  is  the  ideal  gas  constant=8. 0143  J  mol  1  (®K)  »  and  Ta 
is  the  air  temperature  in  °K.  The  quantities  recorded  are  the 
barcaetric  pressure  in  millibars  and  the  air  and  dew-point 
temperatures  in  degrees  Celsius.  Hating  appropriate  unit 
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conversions,  and  substituting  Equation  2.3  in  2.2  gives: 

e=2. 2158*10*  x  (Ta)  x  exp (-5 107. 4/Td)  2.4 

After  obtaining  the  vapour  pressure  Equation  2.4  was  substituted 
in  Equation  2.1  to  qive  the  air  density. 

2.4  Assumptions 

The  nost  sensitive  parameters  in  the  bulk  aerodynamic  flux 
formulae  of  Equations  1.2  are  the  transfer  coefficients.  Por 
example,  the  form  of  the  draq  coefficient  has  undergone 
considerable  debate.  Estimates  vary  from  a  constant  value  (Pond 
et  al.  1971),  to  a  step  function  of  wind  speed  (Rossby  and 
Montgomery  1935),  to  a  linear  function  of  wind  speed  (Smith  and 
Banke,  1975).  Traditionally  aesoscale  and  synoptic  scale 
dynamic  studies  have  employed  the  constant  form.  Recent 

estimates  (eg.  Large,  1979)  at  wind  speeds  greater  than  15 
i/sec  have  indicated  that  the  linear  form  may  be  appropriate. 
Consequently  for  this  study  tvo  drag  coefficient  formulations 
were  used  --  a  constant  fora  of  1.5x10"*  (after  Pond,  1971)  and 
a  linear  fcra  (fcoa  Larqe,  1979): 

Cd  =  1.14x10'*  0<  V  <  10  a/sec  2.5 

Cd  =  (0.  49*0. 065V) X  10'*  V  >  10  m/sec 

where  V  is  the  wind  velocity  averaged  over  one  hour. 

The  Stanton  and  Dalton  numbers  are  less  well  understood 
than  the  drag  coefficient.  Large  lists  instrument  contamination 
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fros  sea  spray  as  a  najor  impediment  in  their  Measurement.  Thus 
a  value  of  1.5*10  ■*  after  Pond  et  al.  (1974)  was  taken. 

Normally  the  transfer  coefficients  are  assused  to  be 
calculated  at  a  reference  level  of  10a  and  under  conditions  of 
neutral  stability.  Neither  of  these  conditions  is  net  by  Most 
of  the  data  used  in  this  study.  The  anenoaeter  level,  for 
exaaple,  at  station  P  is  at  22b.  Over  most  of  the  Northern 
temperate  oceans  the  sea  is  generally  warier  than  the  air  and 
the  sea  surface  saturation  is  generally  greater  than  the 
absolute  huaidity  of  the  surrounding  air  indicating  that 
generally  unstable  conditions  exist. 

Large  (1579)  and  Turner  (  1973)  give  the  corrections  reguired 
to  the  transfer  coefficients  to  account  for  non-noraal 
conditions.  Corrections  applied  in  turbulence  studies  reguire 
the  calculation  of  a  flu*  Richardson  nuaoer  involving  the 
correlation  of  the  vertical  velocity  fluctuations  with  the  air 
temperature  fluctuations  and  the  air  absolute  huaidity 
fluctuations  at  microscale  sampling  rates.  A  bulk  stability 
parameter  could  be  estimated  from  the  data.  Large(1979) 
indicates  that  certain  oceanic  conditions  (eg.  an  air-sea 
teaperature  difference  of  5.0  °C  and  a  wind  speed  of  10.0  a/eec) 
could  induce  about  a  10k  error  in  the  transfer  coefficients  due 
to  stability.  Averaging  the  fluxes  would  lead  to  much  smaller 
errors  and  because  the  study  will  be  concerned  with  ratios  of 
fluxes,  such  errors  will  further  cancel.  further note ,  the 
relationship  between  the  stability  and  the  flux  ratios  at 
Station  A  tor  1949  was  investigated  and  no  correlation  was 
evident.  Thus  stability  considerations  are  lett  as  ainor  errors 
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in  the  data.  Non-standard  height  is  not  corrected  because  it 
vould  involve  a  direct  multiplication  factor  to  the  stresses 
(Pi  ssel  et  al.  1977) .  Since  this  study  will  be  concerned  with 
ratics  of  stress,  the  height  corrections  cancel.  (indeed  any 
constant  formulation  of  the  transfer  coefficients  cancels). 
Thus  it  is  assumed  that  the  transfer  coefficients  defined  in 
Equations  1.3  and  1.4  are  at  10m  reference  height  and  at  neutral 
stabil it  y. 

In  order  to  calculate  the  latent  heat  flux  it  was  necessary 
to  find  the  sea  surface  saturation  and  absolute  air  huaidity 
difference.  The  absolute  air  huaidity  is  given  by  Equation  2.2. 
The  sea  saturation  was  assumed  to  be  at  the  sea  surface 
temperature.  Ts  was  used  in  place  of  Td  to  calculate  the  sea 
surface  humidity.  This  value  was  then  multiplied  by  0.98  to 
account  for  the  depression  of  the  saturation  point  due  to  the 
pre  sence  of  salt. 

Throughout  the  thesis,  the  air-sea  humidity  differences  ana 
net  the  constituent  soa  surface  and  dew  point  temperatures  were 
averaqed.  For  usual  climatic  c  lieu lat ions,  however,  averaged 
humidify  differences  are  seldom  available.  bather,  the  normal 
information  is  *he  averaged  constituent  variaoles.  Consequently 
a  systematic  errci  is  induced  in  most  ca icu la t l ons . 

It  we  consider  a  single  astimate  of  g  to  be  cased  upon  a 
mean  temperature  T  and  a  deviation  from  the  mean  T*  then  (see 
Equation  2.2): 


g  =  C  1  e  x  p  (C2/  ( T  ♦  T  ' )  ) 


2.6 
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Assuming  that  the  mean  temperature  is  much  greater  than  the 
fluctuation,  the  exponential  can  then  be  re-expressed  as  ; 


1 


,  c, 


T' 


which  oecomss: 


V 


2.  7 


A  typical  ratio  of  T 1 2 / T 2  for  the  dew  point  temperatures  was 
found  to  ne  1.  2x10  Allowing  a  fluctuation  of  5  standard 
deviations  and  a  value  of  C2/T=-18.0  (ie.  T=273  °K) ,  the  terms 
above  in  the  exponential  are  correct  to  0(10-*). 

Expanding  Equation  2.7  in  a  Taylor  series  and  keeping  all 
terms  of  0(>10*2)  gives: 


V 
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Averaging  both  sides  of  Equation  2.8  and  recognizing  that 
cl  exp(C2/T)  is  the  unsolute  humidity  determined  from  the 


average  te  it  l  d  t  u  r  1  (q(T))  then  : 


2.9 


If  tne  distribution  is  reasonably  Gaussian  then: 
a .  I • "  - n  odd  =0.0 

c.  TTT=  3x  (T7^)  ^ 

The  terms  in  T '  4/T  4  are  all  0(<10'3)  and  Equation  2.9 
to : 


.  l 


red  uces 
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Taking  (C2/T)  =  -18.0  and  TT7/!  2  =  1.2xl0-4  (a  typical 

value)  we  see  that: 

g  =  q  (T)  (  HO.  02) 

Since  tne  maximum  averaging  period  is  26  days,  an  upper  limit  on 
the  systematic  error  induced  in  averaging  the  humidities  over 
averaging  the  dew  point  temperatures  is  approximately  2.0%. 
1 urther mors,  the  T 1 2/T  2  ratio  for  the  sea  surface  temperatures 
was  about  one  half  that  of  the  dew  point  temperatures.  This 
gives  a  systematic  sea  surface  humidity  bias  of  about  1.0ft. 
bince  the  two  riases  are  suDtracted  in  the  latent  heat  flux, 
this  gives  a  negative  systematic  air-sea  humidity  difference 
mas  of  only  1.0ft.  This  was  considered  negligible  in  comparison 
with  tne  prctatle  sampling  errors  in  the  data. 
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2 . 5  analysis  And  Definitions 

Because  each  ship  contained  a  varying  nuabec  of  years  of 
data,  they  were  reconstructed  into  yearly  blocks  of  2912 
observations  each.  The  2912  figure  was  chosen  because  it  was 
evenly  divisible  by  multiples  of  2,  4,  and  13  allowing  semi¬ 
annual  and  lunar  monthly  divisions  of  the  year.  The  eight  (or 
1b  for  a  leap  year)  observations  at  the  end  of  each  year  were 
omitted  from  calculations.  This  rs  unlikely  to  significantly 
influence  the  final  results. 

From  the  raw  three-hourly  data  set,  the  t  hree -hour  ly 
densities,  absolute  humidities,  stresses  and  heat  fluxes  were 
calculated.  The  sampled  and  calculated  variates  were  then 
averaged  in  sets  of  2,  4,  8,  1b,  12,  5b,  112,  and  224  readings, 
corresponding  to  averaging  over  0.25,  0.5,  1.0,  2.0,  4.0,  7.0, 
14.  C,  and  28.0  days.  Thus  for  a  one  year  saaple,  the  set  of 
variates  averaged  over  224  three-hourly  readings  contains 
2912/224=13  members.  The  length  of  time  over  which  the  three- 
hourly  data  set  was  averaged  will  be  referred  to  as  the 
averaging  period,  will  be  denoted  by  the  capital  letter  L,  and 
will  be  expressed  in  days.  Thus  averaging  224  three-hourly 
measurements  yields  an  averaging  period  of  L=28.l)  days. 

The  fluxes  were  then  recalculated  using  the  averaged 
constituent  data.  The  winds  and  fluxes  calculated  in  this 
manner  will  referred  to  as  the  vector  averaged  or  VA  winds  and 
fluxes.  The  VA  fluxes  were  then  compared  to  the  directly 
averaged  fluxes  from  the  original  three-hourly  time  series.  The 
directly  averaged  fluxes  will  be  referred  to  as  the  three-hourly 
averaged  cr  3H  fluxes. 


An  example  using  the  stress  components  should  clarify  these 


definitions.  Defining  the  positive  x  and  y  directions  as  east 
and  north  respectively,  the  3H  stress  vector  components  averaged 
over  L  days  are: 


(Tx^,Ty^) 


x 

.J-  V 


-  (  \  2-iO*Ld(  (xC  +  u-*-) 
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where  I  =  L  days  x  8  me asur ement  s/da  y ,  u  L  and  vt'  are  the 
constituent  wind  components  from  the  raw  three-hourly  series,  jOc 
are  the  three-hourly  air  densities,  Cd u  are  the  three  hourly 
drag  coefficients,  and  Tx^  and  Ty^  are  the  3H  x  and  y  stress 
components.  Ihe  VA  wind  vectors  over  the  same  averaging  period 
are : 


2.  12 


where  u^  and  are  the  VA  wind  components  averaged  over  L  days. 
A  VA  stress  can  then  be  calculated  from  the  VA  winds  by  : 

2-'3 

^  is  the  density  averaged  over  L  days,  Cd^  is  the  drag 

•  ■ 

coefficient  based  on  the  VA  wind  and  Tx  and  ly,  are  the  VA 

)  a 

stress  component  estimates. 

Tnroughout  the  text,  the  subscript  i  will  refer  to  a  member 
of  the  raw  *hree-hourly  data  set  and  the  subscript  j  will  refer 
to  a  member  of  the  data  set  obtained  from  averaging  the  three- 
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hourly  set  in  L  day  groups.  The  3H  fluxes  will  appear  without 
primes  and  the  VA  fluxes  will  appear  with  primes. 

The  analysis  of  the  differences  between  the  VA  and  3H 
fluxes  requires  three  distinct  steps.  First,  the  raw  VA  and  3H 
variates  must  be  evaluated  statistically  to  determine  the 
inherent  differences  between  the  two  estimates.  Next,  since  the 
3H  variate  is  the  more  accurate  flux  estimate,  transformations 
must  ne  found  for  the  VA  variate  so  that  it  more  closely 
approximates  the  3H  variate.  Once  a  suitable  transformation  has 
been  found,  it  must  be  applied  to  the  VA  variate  and  the 
statistics  re-evaluated.  Thus,  the  analysis  involves  two 
distinct  calculations  —  computing  the  statistical  relationship 
between  the  VA  and  3H  variates  and  computing  the  transformation 
itself. 

Let  Xj  and  X  •  j  be  members  of  the  3H  VA  flux  variates 
respectively  for  an  averaging  period  L.  The  X'j  may  or  may  not 
include  the  transformation  step  and  Xj  and  X'j  define  generally 
tne  x  or  y  components  of  stress  or  the  latent  or  sensible  heat 
fluxes.  The  long-term  means  (T,)^)  are  defined  as: 


(x.K'j 


2.  IV 


wnere  J  is  the  total  number  of  averagings  available  in  the 
record  for  a  given  averaging  period  L.  Similarly  the  long  tern 
variances,  (Tx  and  ,  are  defined  as: 
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The  relationship  between  the  3H  and  VA  variates  is  then 
defined  in  terms  of  four  statistical  quantities  which  will 
collectively  be  called  the  test  functions.  First,  the 
effectiveness  of  determining  the  long  term  mean  flux  is  defined 
by: 


Ax  «  \x-  x\ 
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This  will  be  referred  to  as  the  difference  mean  (or  DM)  •  A 
measure  of  the  conservation  of  the  total  variance,  in  the 
transforma ticn ,  is  given  by: 


atI 
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and  will  be  referred  to  as  the  difference  variance  (DV) .  The 
effectiveness  cf  the  transformation  on  a  pcint-by-point  basis  is 
given  by  the  residual  variance  (EV) ,  defined  as: 


2.  18 


Finally,  the  correlation  coefficient  determines  the  statistical 
dependence  of  one  variate  on  the  other  and  is  defined  by: 
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2.  19 


Assuming  that  the  atmosphere  has  a  climatological  steady 
state,  the  DMs  indicate  whether  the  climatological  mean  of  the 
VA  variate  approaches  that  of  the  2H  variate.  For  spectral 
studies  such  as  Willefcrand  ( 1978)  both  the  total  variance  and  the 
frequency-by-frequency  power  spectral  densities  for  the  two 
variates  should  match.  The  former  condition  is  determined  by 
tne  DVs ;  the  latter  condition  will  be  investigated  in  Chapter  V. 
The  fiVs  determine  the  relationship  between  the  VA  and  3 H 
variates  on  a  pci nt-by-pcin t  basis  (ie.  on  the  time  scales  of 
the  averaging  period  itself) .  Finally,  the  correlation 
coefficients  determine  statistically  the  ability  of  one  variate 
to  infer  the  other.  It  is  quite  conceivable,  for  example,  to 
find  a  scale  factor  which  will  effectively  reduce  the  DMs,  DVs, 
and  aVs,  while  leaving  the  correlation  coefficients  unchanged. 

Ine  DVs,  SVs  and  correlation  coefficients  have  been 
normalized  by  the  3H  variance  for  ease  of  inter-comparison 
between  ships.  This  effectively  removes  the  height  correction 
(discussed  in  Section  2.4)  from  the  statistics  since  it  would 
appear  to  the  same  power  in  the  numerator  and  denominator  in  all 
three  cases.  The  DMs  have  not  been  normalized  because  an 
appropriate  normalizing  factor  would  be  the  mean  3H  flux  value. 
At  certain  locations,  however,  it  approached  0  dPa  for  the 
stress  (Decipascals  are  used  for  units  and  are  numerically 
equivalent  tc  dynes/cm2)  and  0  Watts/mz  for  the  heat  fluxes. 
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The  absolute  value  is  used  in  the  00  calculation  because  we  are 
concerned  that  the  climatological  3H  value  be  estimated 
irrespective  of  sign.  Recent  information  (Large,  personal 
communication,  1980)  indicates  that  the  Stanton  and  Dalton 
nunbers  used  in  this  study  (1.5x10  3)  are  about  20.0%  high. 
Furthermore,  the  non-standard  height  and  the  stability 
assumptions  outlined  in  Section  2.4  aay  induce  a  further  20.0% 
error  in  the  data.  Thus,  when  a  DO  or  a  mean  3H  flux  is  guoted, 
the  values  qiven  aay  be  biased  high. 

A  discussion  of  the  VA  transf or nation  procedures  foras  an 

a 

integral  part  of  the  text  and  will  be  discussed  in  more  detail 
later.  In  all  cases,  however,  the  test  functions  defined  in 
Equations  2. 16  to  2. 19  will  be  used  to  evaluate  the 
eff ecti veness  of  the  transf oraations. 
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CHAPTER  III  MOHJNTUH  fL’JXgS 


The  processed  data  base  outlined  in  Chapter  II  was  analysed 
to  deteraine  the  extent  of  the  expected  unde Lestiaation  of 
stress  when  using  winds  averaged  over  tines  longer  than  one 
hour.  The  reductions  are  studied  both  for  calculations  of 
cliaatological  values  of  long  tern  net  stresses,  (based  on  wind 
data  averaged  over  a  nonth  or  more)  and  for  calculating  actual 
tine  history  values  for  which  data  night  be  averaged  to  about 
one  month.  Climatoloqica 1  stresses  are  used  as  an  input  in 
nunetical  models  attenpting  to  calculate  the  steady  state 
oceanic  circulation.  The  tenporal  variations  are  required  for 
such  studies  as  that  by  Wi llebran d  ( 1978)  for  calculations  of  the 
stress  spectra  usinq  vector  averaged  winds  derived  fron  pressure 
traps  cr  for  models  attempting  to  examine  the  effects  of  time- 
dependent  wind  stress  forcing.  The  bulk,  aerodynamic  formulae 
for  the  calculation  of  the  stresses  are  given  in  Equations 
1.2(a)  and  (b)  . 


3,  1  The  Vector  Averaged  3H  And  V A  Stress 
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The  3H  stress  provides  an  estimate  of  the  'true' 

1  be  examined  first.  The  annual  cycle,  in  terns 
leans  and  standard  deviations  will  be  exanined 
nts  Fer  nonth  givinq  13  lunar  months  per  year. 

Due  to  the  enormous  amount  of  data  analysed,  not 
ults  can  be  presented  in  this  thesis.  All  ships  e 
ter  naxiaum  and  sunnier  minimum  both  in  the  average 
the  stress  variability  of  which  Meatherstation 
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example.  This  pattern  is  indicative  of  a  winter  storm  season 
typical  of  the  northern  temperate  latitudes. 

A  plot  of  the  stress  components  for  Ueatherstation  C 
appears  in  Figures  4(a)  through  (d) .  The  error  bars  are  ♦  1 
standard  deviation  abot’.t  the  lean.  It  is  readily  apparent  that 
the  standard  deviations  nay  be  as  auch  as  one  order  of  aagnitude 
larger  than  the  mean  stress  value  (and  thus  the  variance  to  aean 
squares  ratio  cf  order  100)  in  either  component  for  any  given 
year.  Over  the  whole  record,  the  x  component  mean  value  is 
positive  with  an  overall  average  value  of  aoout  1.0  dPa  (which 
is  numerically  equivalent  to  1.0  dynes  cm  2)  while  the  y 
component  average  value  is  about  0.20  dPa  and  the  x  component 
standard  deviation  is  about  1.25  times  greater  than  the  y 
component  standard  deviation.  The  maximum  difference  in 
standard  deviation  between  years  appears  in  the  x  component 
between  1952  and  1959  where  the  ratio  is  about  2.  Figures  4(c) 
and  (d)  demonstrate  the  annual  cycle  in  variability;  the  ratio 
cf  standard  deviations  between  month  1  (January)  and  month  7 
(July)  is  about  J. 

The  standard  deviation  of  the  monthly  means  was  also 
calculated  and  appears  in  Figure  5  to  show  the  inter-annual 
variability  of  the  monthly  stress  mean.  The  annual  cycle  of  the 
mean  stress  is  evident  particularly  in  the  x  component  which 
varies  from  about  1.5  dPa  in  the  winter  to  about  0.5  dPa  in  the 
summer.  An  annual  cycle  in  the  mean  y  stress  is  not  evident,  k 
comparison  between  the  standard  deviations  of  Figures  5  and  4(c) 
and  (d)  indicates  that  they  have  been  reduced  by  a  factor  of  two 
to  three,  net  the  factor  of  15  one  might  expect  if  the  224 
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Figure  5  (a) 
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Figure  5  (t) 

Figure  5.  The  3H  stress  as  a  function  cf  month  at  Station  C. 
The  solid  line  is  ±  1  standard  deviation  of  the  monthly  means 
using  the  constant  drag  coefficient  and  the  dashed  line  ±  1 

standard  deviation  of  the  monthly  means  using  the  linear  drag 
coefficient.  The  means  are  indicated  by  the  horizontal  bars. 

similar  to  the  time  scales  cf  storms  in  the  mid-latitudes. 

Table  V  gives  the  averages  and  standard  deviations  of  the 
entire  record  for  all  the  stations  examined.  Note  that  Station 
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2ABLE  V 

The  x  and  y  components  of  3H  stress  leans  and  standard 
deviations  foe  all  _ships  and  years  exaamed.  The  standard 
deviation  =  (lk*-Tk2)  ‘/2  where  k  denotes  the  x  or  y  component. 
All  neasurenents  are  in  dPa. 
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N  values  of  the  standard  deviation  are  1/3  to  1/2  of  those  of 
the  ether  weatherships  (or  the  variance  is  1/9  to  1/4  of  the 
ethers) . 

In  Figure  6,  the  3H  stress  components  are  plotted  versus 
the  VA  stress  coaponents  for  the  constant  drag  coefficient  and 
L=28.0  days  for  Station  A.  The  diagonal  solid  line  indicates 
the  equivalence  line  between  the  38  and  VA  stress.  Bather  than 
scattering  around  this  line,  the  absolute  value  of  the  VA  stress 
is  consistently  less  than  that  of  the  3H  stress  although  they 
are  fairly  well  correlated. 

Figure  7,  shows  estiaates  of  the  four  test  functions 
defined  in  Eguations  2.16  -  2.19  for  the  x  coaponent,  constant 
drag  coefficient  when  no  correction  to  the  VA  stress  has  been 


r  COUPON:. \ 


J 


°  \  /  **• 


2  0  6  0 
Vfl  STRESS  IN  DPfl 


2  0  6  0 
VR  STRESS  IN  DPR 


Figure  6.  3H  versus  VA  stress  for  the  constant  drag  coefficient 
at  Station  A,  L=28.  0  days  with  no  corrections  applied. 


applied.  The  actual  values  for  both  components  and  drag 
coefficients  are  given  in  Appendix  B  .  Symbols  that  will  be 
used  in  the  figures  throughout  the  text  to  indicate  ship  are 
given  in  Table  VI. 

In  Figure  7(a),  the  difference  means  are  small  for  short 
averaging  periods  indicating  that  the  VA  stress  closely 
approximates  the  long  term  3H  mean.  As  the  averaging  period 
increases,  the  VA  estimates  become  increasingly  inaccurate  and 
the  errors  appear  to  increase  approximately  exponentially  with 
averaging  period.  There  is  a  distinct  ordering  in  the  DUS  by 
ship  through  all  averaging  periods.  A  comparison  of  Appendix  B 
and  Table  V  indicates  that  the  ordering  may  be  dependent  upon 
the  mean  3H  stress  component  value  averaged  over  the  entire 
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Figure  7(a)  DMs  in  ifa 


Figure  7(d) 


Figure  7(c)  Figure  7(d) 

Figure  7.  Example  of  the  composite  test  functions  for  all 
ships,  x  component,  constant  draq  coefficient  with  no  correction 
applied . 

record.  Ihose  ships  with  larger  initial  3ri  values  appear  to 
have  large  difference  mean  stresses. 

ln=  difference  variances  appear  in  Figur®  7(d).  For 
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TABLE  VI 

Syabcls  used  in  the  figures  throughout  the  text. 
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Note:  When  it  is  indicated  that  Station  N  is  excluded  from  the 
fiquLe,  Station  E  is  indicated  by 

reference,  the  1H  stress  variances  for  all  averaging  periods 
appear  in  Table  VII.  Again  the  ?A  variance  closely  approximates 
the  JH  variance  for  small  averaging  periods.  By  L=1.0  days, 
however,  the  VA  variance  underest i nates  the  JH  variance  by  about 
.0*.  At  1  =  28  days,  the  variance  reduction  is  60  to  80%.  The 
ship-tc-ship  DVs  are  more  tightly  grouped  than  the  difference 
tea  nr.  The  relitivt  position  of  each  ship  appears  to  shift  as 
the  averaging  period  increases. 

Th»  residual  variances  appear  in  Pigure  7(c).  At  L-0.  2 5 
days,  the  residuals  are  less  than  1.0%  of  the  total  3H  variance. 
Py  1=4. c  days  this  i lses  to  a  lean  of  aoout  10%  and  by  20  days 
the  value  is  between  20  and  49J  of  the  total  3d  variance. 

overall,  the  correlation  coefficients  shown  in  Figure  7(d) 
are  quire  high.  At  L=c. 2S  days,  all  values  are  greater  than 
0.99  while  at  L=28.0  days  all  values  lie  between  0.93  and  0.97. 

Alrhouqh  only  the  i  coiponent,  constant  drag  coefficient 
appears  here,  the  lanqe  ot  values  tor  the  y  component  constant 
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tabu  111 

The  stress  variances  for  both  drag  coefficients  as  a  function  of 
averaging  period.  The  columns  refer  to  the  averaging  length  in 
days  and  the  variances  are  in  (dPa) *. 
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TABLE  VII  (con  tinued) 
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dray  coefficient  are  similar  to  those  quoted  above  with 
differences  occurring  only  in  the  detail  of  the  scatter.  The 
difference  means,  difference  variances  and  residual  variances 
are  greater  for  all  ships  and  averaging  periods  using  the  linear 
drag  coefficient  than  when  using  the  constant  drag  coefficient 
while  the  correlation  coefficients  are  consistently  smaller. 


3.2  Previous  Studies 


Fissel(1975)  analysed  10  years  of  data  from  Station  P.  He 
divided  his  data  into  five  two-year  olochs  of  5832  data  points 
each  and  quantified  the  reduction  in  terms  of  two  functions 
which  I  define  as  : 
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wnere  T'x  and  I '  y  are  the  x  and  y  component  VA  estimates  of  the 
two-year  mean  stresses  for  an  averaging  period  of  L  days  and  Tx 
and  ~if7  are  the  two-year  mean  3H  stress  components.  One  S  (L)  and 
A  (L)  value  was  calculated  for  each  5832  point  block,  and  the 
sample  means  and  standard  deviations  were  calculated  over  the 
five  blocks. 

The  Fissal  results  appear  in  Figure  8.  (Fissel  used  the 
linear  Cd  form  cf  Deacon  and  Webb(1S£2)  rather  than  the  one  used 
here).  The  95%  confidence  zones  (equal  tot  2  standard 
deviations  of  the  two-year  means)  were  typically  of  the  order 
0.1  times  the  mean  value  cf  S(L)  and  2°  for  A (L) . 

S  (L)  determines  the  ratio  of  the  climatological  average  VA 
stress  magnitude  to  the  average  3H  stress  magnitude  and  A  (L) 
determines  the  average  difference  between  the  3H  and  VA  stress 
directions.  Exact  long-term  mean  3H  stress  from  VA  stresses 
within  the  two-year  averaging  block  should  result  from  the 
application  of  S(L)  and  A  (L)  to  the  VA  stresses.  The  test 
functions  were  calculated  by  applying  S  (L)  and  A(L)  to  the  data 
and  the  DMs  were  substantially  improved  but  the  other  test 
quantities  remained  nearly  identical  to  those  in  Figure  1  where 
no  correction  was  applied. 

Tnis  technique  was  attempted  on  the  North  Atlantic  ship 
data.  In  the  present  study,  it  was  found  that  one-yearly  blocks 
of  2912  data  points  simplified  the  analysis  when  using  data  of  a 
varyiny  number  of  years.  For  most  weatherstations,  the  Station 
P  results  were  well  replicated  as  shown  in  Figures  9(a)  and  (fc) 
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for  Station  C.  The  S (L)  values  at  Weatherstation  A  (Figures  9 
(c)  and  (d) ) ,  however,  differ  markedly  froa  the  other  ships. 
Whereas  Stations  P  and  C,  S(L)  values  are  consistently  less  than 
1.0,  Station  A  has  several  values  significantly  greater  than 
1.0.  The  standard  deviations  for  averaging  periods  froa  0.5  to 
28.0  days  are  as  auch  as  100  tines  larger  at  A  than  C  and  are 
large  compared  to  the  nean  S (L)  values.  Station  A  has  direction 
shifts  ( A  (L) )  of  up  to  30°  for  averaging  periods  greater  than  14 
days.  For  a  large  value  to  arise  in  S  (L) ,  the  denominator  of 
Equation  3. 1  nay  tend  to  0  while  the  nuaerator  reaains  finite. 
One  would  expect  high  values  of  S  (L)  to  occur  for  low  values  of 
the  average  3H  stress  magnitude.  The  yearly  values  of  S  (L)  at 
Station  A  are  listed  in  Table  VIII  along  with  the  yearly 
averaged  3h  stress  Magnitudes. 

From  Table  VIII  it  is  clear  that  generally  low  yearly 
averaged  stress  values  lead  to  S  ( L.)  >  1  (also  true  for  the  linear 
drag  coef f icient) .  This  is  not,  however,  a  sufficient 
condition.  The  1955  stress  naqnitude  is  0.3  dPa  while 
S  (0. 25)  =  1 . 0005 ;  in  1952  the  stress  is  0.107  dPa  but  S(0.25)<1.0. 
The  highest  stress  value  for  which  S  (L)  is  greater  than  1.0  is 
at  Weatherstation  B  in  1965  where  the  stress  was  0.372  dPa.  In 
general,  hcwever,  saall  vector  averaged  3H  stress  magnitudes 
lead  to  large  values  of  S  (L)  .  Furtheraore  the  lower  the  3H 
stress,  the  nore  averaging  periods  are  affected.  All  values  of 
S(L)  greater  than  2.0  have  vector  averaged  3H  stress  magnitudes 
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Figure  9.  S  (L)  and  A  (L)  (in  Degrees)  as  functions  of  averaging 
period  for  Weatherships  A  and  C.  The  mean  for  the  constant  drag 
coefficient  is  indicated  by  A  with  ±  1  standard  deviation  of 
che  yearly  means  by  the  narrow  horizontal  lines.  The  mean  for 
the  linear  drag  coefficient  is  indicated  by  >c  and  ±  1  standard 
deviation  of  the  yearly  means  by  the  wide  horizontal  lines. 
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TABLE  VIII 

The  S  (L)  readings  for  all  years  of  data  at  Weatherstation  A 
using  the  constant  drag  coefficient  .  The  3H  vector  averaged 
stress  aaqnitude  is  in  dPa. 
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3*3  Individual  Ratios 

In  attempting  to  quantify  the  reduction  in  stress  due  to 
vector  averaging  the  wind,  the  functions  S (L)  and  A  (L)  have 
several  disadvantages.  First-,  in  regions  of  small  average 
stress,  the  functions  may  assume  radically  different  values  from 
one  year  to  the  next,  indicating  that  they  may  not  be  generally 
applicable  to  years  other  than  these  over  which  the  initial 
averaging  was  dene.  Second,  the  S(L)  and  A  (L)  functions  do  not 
improve  the  PVs  and  correlation  coefficients.  For  calculations 
of  actual  time  variations  it  is  necessary  to  estimate  the 
accuracy  of  the  individual  stress  values.  Consequently  two 
functions  B1(l)  and  Bj(L)  were  defined  as: 


3.3 


wnere  T»Xj-  and  T'y^  are  the  sets  of  VA  x  and  y  component  stress 
estimates  and  Tx^  and  Tyj . are  the  sets  of  3H  x  and  y  component 
stress  estimates  for  an  averaging  period  of  L  days  following  the 
notation  given  in  Section  2.5. 


The  Rj(L)  for  any  given  averaging  period  can  be  examined 
mathematically  to  gain  some  insight  into  its  expected  behaviour. 
Eguation  3.3  for  the  constant  drag  coefficient  can  be  rewritten 


witn  the  approximation  that  the  air  density  and  the  drag 
coefficient  are  constant.  Applying  the  Schwartz  ineguality  (see 
Pfeiffer,  1965,  pp.  223)  to  the  denominator  and  rearranging 


terms,  a  lower  bound  for  the  function  can  be  established: 


_1 _ 

-  St 


Where:  and  CT>^'  are  the  x  and  y  wind  velocity  component 
variances  in  averaging  period  j.  Thus  the  lower  bound  for  any 
particular  averaging  period  increases  for  higher  vector  averaged 


mean  winds  and  lower  variances 
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Rj(L)  has  no  upper  bound.  If  one  takes  the  set  u  =  v  * ( 
3, 4, -5  }  and  substitutes  into  Equation  3.5,  then  the  numerator 
is  finite  and  the  denoninator  is  0.0,  or  Rj(L)  =  «0  .  The 
probability  of  such  an  occurrence  is  reaote.  Whenever  the  3H 
stress  averaged  exactly  to  0.0,  leading  to  fij(L)=QD  ,  the 
particular  averaging  values  were  eliminated  froe  calculations. 
After  imposing  this  criterion,  in  over  750,000  cases  (all  ships, 
averaging  periods  and  drag  coefficients)  examined,  not  one 
instance  of  R j (L)  greater  than  1.0  appeared.  Furthermore,  in 
all  cases  where  the  3H  stress  was  0,  the  mean  VA  stress  was  also 
0. 

The  smoothed  wind  data  as  encountered  in  surface  pressure 
maps,  hides  the  internal  variance  through  the  averaging  process, 
hut  the  square  of  the  vector  averaged  wind  is  not  only  readily 
calculable  but  essential  in  determining  the  value  of  VA 
fluxes.  Although  the  value  of  any  single  Bj(L)  and  the  lower 
bound  nay  net  coincide,  an  average  value  (8  (L) )  may  be  wind 
speed  dependent.  Consequently  the  effect  of  grouping  Hj(L) 
value  according  to  the  vector  averaged  wind  speed  was 
investigated. 

F j f L)  and  Bj(L)  were  sorted  into  twenty  equally  spaced 
divisions  according  to  value.  For  each  Rj(L)  and  B j (L)  reading 
the  vector  averaged  wind  speed  (u  v^ 2)  l/2  was  determined  and 
the  functions  were  further  sorted  into  thirteen  divisions 
according  to  the  Beaufort  number  of  the  j  value  of  the  vector 
averaqed  wind  speed.  The  Beaufort  wind  speeds  according  to  the 
Narine  Climatic  Atlas,  U.S.  Office  of  Naval  Operations  are 
given  in  Table  IX. 
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TABLE  IX 

BEAUFORT  WIND  SPEED  INTERVALS 
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Figure  10  is  the  histograa  of  Rj(L)  and  Figure  11  the 
histogram  for  Ej(L)  for  the  constant  and  linear  drag 
coefficients  respectively  at  Weatherstation  C.  Beaufort  nuabers 
2  (1.C  a/sec),  6  (9.4  a/sec),  and  10  (22.65  a/sec)  for  an 
averaging  period  of  L=0.25  days  were  selected  as  being 
representative.  The  standard  deviation  was  calculated  froa  all 
the  values  of  Rj(L)  and  B j  (L)  in  each  wind  speed-averaging 
period  category  and  will  be  referred  to  as  the  saaple  standard 
deviation  and  its  square  will  be  referred  to  as  the  saaple 
var  iance. 

As  seen  in  Pigure  10  at  1.0  i/sec,  Hj(L)  has  a  distinct 
■ode  between  0.45  and  0.50  with  an  average  value  of  0.530.  The 
rest  of  the  distribution  is  nearly  flat.  At  9.40  a /sec  the 
distribution  of  8j(L)  becomes  markedly  skewed  towards  1.0.  The 
standard  deviation  has  also  significantly  decreased.  At  22.65 
■/sec,  the  skewness  becomes  even  lore  pronounced  to  1.0  with  93* 


'  f.  *-■ 


DEGREES 


4  9 

of  the  distribution  in  the  0.9  5  to  1.0  slot.  Because  of  me 
increasing  concentration  of  values  to  the  nigh  end  of  the 
distxin ution ,  the  standard  deviation  continues  to  decrease  with 
averaging  speed.  The  trend  to  higher  R  (L)  values  and  lower 
standard  deviations  for  increasing  wind  speed  is  generally  true 
for  all  averaging  periods. 

Bj  (L)  has  a  distinctly  spiked  distribution  at  low  wind 
speeds.  The  histogram  for  1.0  m /sec  has  £4%  of  the  total  number 
of  3t>3  points  in  the  0  to  1.0°  direction  shift  slot.  There  are, 
however,  an  additional  59  points  from  -10°  to  -180°  and  an 
additional  68  points  from  10°  to  180°  accounting  for  the  large 
18.021°  standard  deviation.  As  the  wind  speed  increases  to  9.40 
m/sec,  the  mean  remains  very  close  to  0.0°.  The  number  of 
exterior  points  decreases  causing  the  sample  standard  deviation 
to  decrease  so  that  at  22.65  m/sec  the  distrinution  appears 
approximately  Gaussian. 

The  joint  histogram  for  Rj(L)  and  B j  (L)  was  calculated  for 
all  snips.  Ihe  arithmetic  means,  1T(T)  and  B (L) ,  were  then 
calculated  for  each  Beaufort  category  and  averaging  period. 
Thus  each  ship  was  assigned  R  (L)  and  1T(L)  values  that  were 
functions  of  windspeed  and  averaging  period. 

Tne  mean  E  (L)  values  for  all  the  ships  appear  in  Figure  12 
while  the  means  E  (L)  values  appear  in  Figure  14.  Two  standard 
deviations  are  also  shown.  Particularly  fcr  larger  VA  wind 
speeds  and  longer  averaging  periods,  R  (L)  and  B  (L)  may  be 
calculated  from  differing  numbers  of  measurements.  First, 
weights  proportional  to  the  number  of  estimates  within  each 
Beaufor t-averaging  period  category  were  applied  to  the  ship 
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vilijf.s  ct  b  (L)  and  to  ar  r  i  vp  at  an  estimate  of  the 

-  tn.iai  i  i-viat  :c!.s  jl  trie  geographic  means.  This  is  the 
:ar  i  v  iat  10:.  at  the  .nears  snown  in  the  figures  and  it  is 

i ..  1  ij  f  U  lushed  lines.  Next  similar  weights  were 

■M  i  1  r  -  to  tae  sample  var i  mctrs  to  arrive  at  an  average  sample 
vaiiu.'.e-1.  I:. is  indicative  of  tne  average  sample  variability 
a.iu-  a  a  y  l •  1 1  -=  f<  (L )  *- sr  i  m are  and  is  indicated  by  the  solid 
1-:.*  -n  Frjures  12  and  14.  Finally,  it  will  be  shown  that 
atut-oi.  N  lifters  markedly  from  other  locations  and  its  values 
ar>.  •  x<..l  id‘  1  entirely  f ran  star,  lard  deviation  estimates. 


hr  Hr;  12  (  i) 

snows  K  (0.  Lb 

)  (tne 

0.25 

indicates  L=0. 25  days) 

i  ur.ct  ion  of 

wind  i pee  ]. 

Tne  v  j 

lue  of 

R  (0.25)  approaches  0 

vtry  small 

wind  speeds. 

F  or 

larger 

vector  averaged  winds 

n(j.it)  rapidly  increases.  At  1  m/sec  R  ( 0 . 2  5)  is  0.90.  As  the 
wij.J  speed  increases,  the  values  approach  1.0.  Both  the 
standard  deviations  of  the  sample  and  the  standard  deviations  of 
the  snips*  means  have  maxima  at  low  wind  speeds  and  then 
decrease  gradually  as  the  wind  speed  increases.  For  all  wind 
speeds,  the  average  error  in  the  sample  is  much  larger  than  the 
error  in  the  means  indicating  that  the  value  cf  R  (0 .25)  may  be 
independent  of  geographical  location. 

Ihe  results  can  be  explained  in  terms  of  Equation  3.6.  An 
averaging  period  of  0.25  days  contains  only  two  three-hourly 
measurements.  Lower  vector  averaged  winds  result  from 
simultaneous  cancellation  of  the  x  and  y  wind  components.  Thus 
tne  variances  will  be  much  greater  than  the  means  of  the 
components  resulting  in  a  lower  bound  that  approaches  0.  For 
higher  wind  speeds  one  expects  a  strong  steady  wind  to  blow  for 
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Figure  12  (c) 

Figure  12.  Piets  cf  R  (L)  for  representative  averaging  periods 
and  wind  speeds  of  all  ships  using  the  constant  drag 
coefficient.  Y  denotes  Weatherstation  N .  The  solid  line  is 
a  weighted  mean  sample  standard  deviation  and  the  dashed  line  is 
the  standard  deviation  of  the  ships’  means.  Since  N  differs 
from  the  other  ships,  its  values  were  not  included  in  the 
calculation  cf  the  standard  deviations. 
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a  relatively  long  length  of  time  in  which  case  the  variance  of 
the  two  components  is  small  compared  to  the  vector  averaged  wind 
and  the  lower  hound  in  Equation  3.6  approaches  i.  R(4) 
approaches  0  as  the  wind  speed  goes  to  C  m/sec.  The  value  of 
i.  (h)  increases  with  wind  speed.  The  scatter  in  measurements  for 
uoth  the  mean  and  sample  values  is  larger  than  for  R(0.25).  At 
the  largest  wind  speed,  the  scatter  in  the  means  appears  to  be 
mere  significant  than  the  scatter  in  the  individual  estimates. 
Die  roar  values  at  21  m/sec,  however,  comprise  less  than  1"  of 
the  total  number  cf  points  measured  and  this  value  cannot  be 
considered  statistically  reliable.  From  0  to  14  m/sec, 
Weatnersnip  N  has  systematically  higher  values  than  the  others. 

Figure  12(c)  exhibits  further  evidence  of  the  higher  values 
of  ML)  for  Weathership  N.  Here,  E  (L)  for  varying  averaging 
periods  is  plotted  for  a  constant  wind  speed  of  6.75  m/sec.  The 
Station  N  R  (L)  values  are  more  than  two  standard  deviations  of 
tne  mean  larger  than  any  ether  ship  indicating  that  its  values 
are  significantly  larger  than  the  ether  stations.  The  scatter 
in  the  individual  measurements  reaches  a  maximum  for  an 
averaging  period  of  two  days  and  then  decreases  slowly  with 
increasing  averaging  period.  The  scatter  in  the  ships'  means 
increases  rapidly  tc  two  days  and  then  maintains  a  relatively 
constant  value. 

Ine  anomaly  in  R  (L)  at  Station  N  is  consistent  both  with 
Equation  3.6  and  the  findings  of  other  researchers, 
idalkus  (  196  2)  notes  that  Station  N  is  at  the  northern  extremity 
or  tne  Northern  Pacific  Trade  Hind  region  characterized  in  both 
direction  an!  magnitude  by  moderate  and  steady  winds.  This  is 
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confirmed  in  Tatle  V  where  the  total  3H  variance  is  about  1/10 
times  that  of  the  more  northerly  stations.  Thus  Equation  3.6 
predicts  that  one  should  get  systematically  larger  values  of 
h j  (L)  and  hence  R  (I) .  This  is  confirmed  by  Malkus  where 
substantially  smaller  discrepancies  between  the  3H  and  VA  stress 
over  a  single  monthly  averaging  were  found  in  the  Caribbean  Sea 
at  19°  N,  66°  H  than  were  found  at  Station  C.  Thus  the 
corrections  required  for  the  tropics  and  subtropics  may  oe  much 
less  than  those  required  for  temperate  latitudes. 

Figure  13  shews  the  geographic  (ie.  average  of  all  the 
stations)  average  values  of  R  (L)  for  the  constant  drag 
coefficient,  excluding  N ,  as  a  function  of  wind  speed.  The 
weighting  outlined  above  was  employed  in  the  calculation  of 
these  geographically  averaged  values.  All  the  curves  appear  to 
be  regular  with  the  exception  of  the  highest  wind  speeds.  This 
may  be  a  result  of  there  being  too  few  averagings  at  the  higher 
wind  speeds  to  yield  a  representative  value  of  R(L). 

The  scatter  and  cehaviour  of  E  (1)  as  a  function  of 
averaging  period  are  indicated  in  Figure  14.  There  is  a  high 
aegree  of  scatter  at  a  wind  speed  of  0.2  m/sec  for  both  the 
sample  and  between  ships.  There  does  not  appear  to  be  a 
systematic  difference  between  ships.  At  6.75  m/sec,  the  sample 
standard  deviation  is  reduced  from  35°  to  about  15°  and  the 
variability  between  ships  is  nearly  zero.  Figure  14(c)  clearly 
indicates  tnat  non-zero  values  of  E  (L)  occur  only  for  low  wind 
speeds.  For  winds  greater  than  about  3  m/sec,  B(L)  averages  to 
exactly  zero.  It  is  net  clear  whether  the  large  values  of  3 ( L ) 
for  low  wind  speeds  are  statistically  significant  since  the  low 
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Figure  13.  Geographically  averaged  B  (L)  values  for  all  wind 
speeds  and  averaging  periods.  lie  numbers  indicate  lines  of 
equal  averaging  periods  (in  days).  Station  N  has  been  excluded 
from  the  averaging. 


wind  speeds  account  for  only  a  very  small  percentage  of  the 
total  number  of  wind  measurements  made. 


3.4  Multiple  Regression  Analysis 

From  the  previous  discussion  it  is  evident  that  the 
reduction  in  stress  due  to  vector  averaging  of  winds  varies  with 
the  averaging  period  and  the  measured  VA  winds.  For  wind  speeds 
greater  than  3.0  m/sec  the  average  shift  in  direction  is  0° 
indicating  that,  on  average,  there  may  not  be  significant 


differences  between  the  transformations  required  for  the  x  and  y 
components.  This  implies  that  a  single  transformation  based  on 
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the  stress  magnitude  is  feasible.  A  single  set  of  correction 
factors  equally  applicable  to  both  components  reduces  tne 
required  parameters  by  a  factor  of  2. 

Tne  object  of  any  transformation  is  to  reduce  the 
difrerences  means,  difference  variances,  and  residual  variances 
to  predict  accurately  the  long-term  and  point -by- point  3H  stress 
from  the  VA  stress.  Specifically,  if  the  point-by-point  3H 
stress  (as  quantified  by  the  residual  variance)  can  be 
accurately  determined  then  the  difference  variances  and  means 
snouid  also  be  reduced.  The  residual  variance  can  be  optimally 
reduced  by  a  least  squares  fit. 

To  this  end  the  3H  and  VA  stress  magnitudes  were  calculated 
and  sorted  according  to  the  Beaufort  interval  as  previously 
outlined.  For  each  averaging  period,  the  following  quantity: 


was  minimized  by  : 

- '  -  o 

c)  3.  8 


giving  : 


'tte 


T1 


£  T^7- 

J  * 

where  T  and  T’v^  are  the  3H  and 
the  stress  magnitude  residual 
period,  l is  a  multiplicative 


3.  9 


VA  stress  mag ni tudes  ,  £><t  is 

variance  for  the  1  averaging 
factor  which  minimizes  the 
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residuals  in  Beaufort  category  k  and  the  1  averaging  period, 
N  is  the  total  number  of  averagings  available  in  the  record  and 
Jtt  is  the  number  of  averagings  within  Beaufort  category  k. 

is  the  slope  of  a  stress  magnitude  regression  line 
passing  through  the  origin.  This  is  not  a  bad  assumption  in 
light  cf  Figure  6,  where  the  scatter  of  the  x  and  y  stress 
components  neccmes  constricted  at  the  origin.  Appendix  A  gives 
an  outline  of  seme  of  the  relationships  between  the  difference 
variances,  residual  variances,  and  correlation  coefficients  for 
the  regression  analysis.  Notice  that  Equation  3.7  unlike 
Equation  2.18  excludes  the  mean  terms.  They  were  excluded 
because  as  demonstrated  in  Figure  4  and  Table  V,  the  stress 
component  fluctuations  dominate  the  mean.  Appendix  A. 1  includes 
the  influence  of  the  linear  term  on  the  residuals.  This  form  of 
the  regression  will  be  applied  presently  tc  the  heat  fluxes.  In 
fact  a  regression  cf  the  type  outlined  in  Appendix  A.  1, 
including  the  mean  stress  magnitude  was  tried  on  the  stress 
magnitudes.  Up  tc  Beaufort  interval  3,  seme  of  the  £*«.  values 
were  found  tc  be  negative  in  contradiction  with  Eguation  3.3. 
Furthermore,  the  component  residual  variances  were  significantly 
larger  than  when  the  means  were  excluded. 

An  estimate  of  95X  confidence  zone  of  the  slope  was  made 


3.  10 


where  tk,  oiS"  is  the  Student  t  statistic  with  J  -1  degrees  of 
freedom  it  a  0.S5  level  of  significance.  This  is  not  a  precise 
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estimate  of  the  error  because  as  shewn  in  the  previous  section, 
the  distribution  cf  Fj(L)  is  highly  skewed  around  the  mean  fer 
most  averaging  periods  and  wind  speeds  and  consequently  the 
statistics  are  not  Gaussian. 

3.5  Variance  Underestimation 

The  values  were  calculated  for  each  ship  and  the  four 
test  functions  outlined  in  Section  2.6  were  then  found  by 
applying  to  the  VA  stress  compon  snts. 

The  difference  variances  using  the  constant  drag 
coefficient  are  shown  in  Figure  15.  Values  for  all  ships, 
components  and  drag  coefficients  are  greatly  reduced  from  the 
values  shown  in  Figure  7.  All  values  are  positive  definite 
which  implies  that  the  long  term  3H  variance  is  systematically 
larger  than  the  long  term  corrected  VA  variance.  Although  only 
tne  constant  drag  coefficient  results  are  shown,  every  case 
examined  showed  a  positive  difference  variance. 

This  is  an  artifact  cf  the  regression  analysis.  Appendix 
A. 1  shows  that  when  the  optimal  is  applied,  the  difference 
variance  and  the  residual  variances  must  be  identical.  Note 
that  even  though  Appendix  A. 1  applies  to  including  the  long-term 
difference  means  in  the  calculation  of  ,  identical  results 
can  be  obtained  by  setting  the  DMs  to  0.  Since  the  residual 
variance  must  be  positive  definite,  so  is  the  difference 
variance.  For  this  result  to  be  exactly  applicable  in  the 
stress  components,  the  x  and  y  components  must  be  regressed 
independently  and  a  separate  set  cf  transformations  found  for 
each  component.  Here  the  stress  magnitude  was  regressed  and 
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Figure  15.  Systematic  variance  underestimation  for  ail  ships, 
constant  drag  coefficient  using  the  initial  estimates  of 
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applied  identically  to  each  component.  Thus  the  applied 

transformation  may  differ  from  the  optimal  component 

transformation.  The  fact  that  the  DV  corrections  in  this 
section  are  positive  definite  indicates  that  the  calculated 
Zvie  based  on  the  stress  magnitudes  approximates  the  optimal 
component  corrections  reasonably  well.  Finally,  this  suggests 
that  a  trade-cff  occurs  in  which  one  can  either  reduce  the  RVs 
to  a  minimum  value  and  accept  a  systematic  error  in  the  long¬ 
term  variance  or  minimize  the  DVs  and  accept  an  increase  in  the 
res idua Is . 

If  the  difference  variances  for  each  component  were 
identical  then  this  bias  could  be  made  exactly  0  by  multiplying 


^  /<r> 


Tv  I  CT\ 


/ 


by  either  x  /fy  or  Tv  j  where  (  ,  <T\y  ) 

and  (  (T* 1  ,  1  )  are  the  component  3H  and  VA  standard 
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deviations  with  the  transformation  applied  to  the  VA  variate. 
This  is  aot  exactly  the  case.  Consequently  a  new  E ^  term  was 
round  oy  : 


V  /V 
£  1  4/ 


3.  1  1 


Tne  has  teen  multiplied  by  the  mean  correction  required  for 

ooth  tne  x  and  y  VA  stress  components. 

Tne  D  Vs  fcr  the  corrected  tran sf or  mat icn  appear  in  Appendix 
0,  and  the  specific  constant  drag  coefficient  results  will  be 
shown  later.  The  variance  differences  are  now  randomly 
scattered  positive  and  negative.  Appendix  D  shows  that  for  all 
ships,  components,  drag  coefficients  and  averaging  periods,  a 
positive  difference  variance  in  the  x  component  implies  that  the 
y  component  DV  is,  within  several  tenths  of  a  percent,  the  exact 
negative. 

Since  the  objective  has  been  to  minimize  the  RV,  the 
correction  applied  to  £  to  remove  the  systematic  bias  from  the 
variance  differences  perturbs  the  residual  variance  so  that  it 
is  no  longer  a  minimum.  The  amount  cf  this  perturbation  can  be 
calculated  under  two  different  sets  of  assumptions  as  outlined 
in  Appendices  A. 2  and  A.  3.  If  we  define: 

Cl  te')  =  £*.*/  L  3. 12 

and  assume  that  (Tx/uV--  ^  t  then  it  is  demonstrated  in 


Appendix  A. 3  that: 


61 


T_  £_  {(£SfeO/C2  0veVl-KO  ±€']\ 

<LT 


3.  13 


wnere  S*-t  is  the  rv  after  (1+£)  has  been  applied.  The 
relative  increase  in  the  RV  can  be  determined  ir.  terms  of  the 
applied  correction  ana  the  correlation  coefficient.  Figure  16 
shows  this  function  for  the  ranges  of  <=  and^T  associated  with 
this  problem.  Excluding  Station  N ,  Figure  15  indicates  that  the 
largest  DV  values  are  about  0.C8  which  gives  a  typical  value  of 
0.04  for  &-  .  A  lew  estimate  of  the  correlation  coefficient 
is  0.95.  When  Equation  3.  13  is  then  evaluated,  the  relative 
increase  in  RV  is  2.4$.  Allowing  for  a  residual  variance  as 
high  as  50$,  the  aosolute  change  in  FV  is  aoout  1%.  (The 
initial  residual  variances  are,  in  fact,  significantly  less  than 
10>u).  The  actual  change  in  RV  was  calculated  and  in  all  cases 
was  substantially  less  than  1.0$. 

henceforth,  rhe  primes  will  be  dropped  and  it  will  be 
understood  that  includes  the  DV  correction. 

3.6  Wind  Dependent  Correction  Factors 

Tue  correction  factors  for  all  ships,  averaging  periods 
and  wind  speeds  are  given  in  Appendix  C  -  Part  I.  The 
corresponding  error  estimates  of  the  95$  confidence  intervals  on 
these  values  as  estimated  by  the  Student  't'  test  are  given  in 
Appendix  C  -  Part  II.  Figure  17  shows  the  inverse  of  correction 
factors  (  l/£*e)  as  functions  of  wind  speed  and  averaging  period. 
The  inverse  is  plotted  for  ease  of  comparison  with  Figure  12. 
Also,  the  inverse  ratio  is  hounded  between  0  and  1.0.  Since  £b<, 
increases  as  the  wind  speed  approaches  0  m/sec,  the  inverse 
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3.7  Accuracy  Cf  Iransfcriaations 

Tn^  difference  means,  difference  variances,  residual 
variances  and  correlation  coefficients  were  determined  by 
applying  and  the  results  appear  in  Appendix  D.  These 
quantities  will  be  examined  in  light  of  their  actual  values  and 
the  arithmetic  difference  between  the  unimproved  values  of 
Appendix  E  and  the  improved  values  cf  Appendix  D. 

As  an  example,  the  difference  means  for  the  x  component, 
constant  drag  coefficient  appear  in  Figure  18(a).  A  comparison 
with  Figure  7  illustrates  several  interesting  features.  First, 
the  difference  means  have  teen  reduced  by  a  factor  of  about  10. 
Second,  there  appears  to  be  an  inverse  relation  in  the  amount  of 
reduction  achieved.  In  Figure  7  Stations  A  and  N ,  with  no 
correction  applied,  have  the  smallest  difference  means.  In 
figure  18,  after  application  of  ,  they  now  have  the  largest 
difference  mean  values.  Conversely,  Stations  C  and  D  have  the 
highest  difference  means  with  no  correction,  but  are  less  than 
0.02  dFa  for  all  averaging  periods  after  correction.  The 
improvements  appear  in  Figure  19(a)  for  the  x  component, 
constant  drag  coefficient.  In  all  cases  there  was  some 
improvement.  An  inspection  of  Appendices  8  and  D  indicates  that 
the  most  dramatic  improvements  occur  where  the  uncorrected 
values  were  initially  large  --  in  particular  for  longer 
averaging  periods  and  the  linear  drag  coefficient. 

The  difference  variance  values  are  now  randomly  scattered 
and  Dcunded  between  about  ±  6 i  using  the  constant  drag 
coefficient  and  about  ±  10X  using  the  linear  drag  coefficient  -- 
again  a  10  fold  improvement  over  the  uncorrected  values.  The 
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results  after  application  of  the  individual  ships'  £ne 
correction  factor  for  the  x  component  constant  drag  coefficient 
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most  accurate  values  occur,  not  unexpectedly,  for  the  lower 
averaging  periods.  The  differences  between  Figures  18  and  7  for 
each  averaging  period  appear  in  Figure  19(b).  Again  all  changes 
are  positive,  indicating  that  the  corrections  have  improved  the 
long-term  variance.  Indeed  for  L  =  28  days  there  are  improvements 
up  to  85%  in  the  3tf  variance. 

Figure  18(c)  shows  the  final  calculated  residual  variances 
for  the  x  component  constant  drag  coefficient.  Note  that  the  RV 
axis  scale  has  been  increased  by  a  factor  of  5.0  from  Figure 
7(c).  With  the  exception  cf  Station  N  all  values  are  less  then 
6 Si.  Figure  IS  ( c )  indicates  the  improvements  achieved.  For  low 
averaging  periods,  the  differences  are  quite  small.  By  L=4 
days,  however,  an  improvement  of  about  5%  has  been  achieved. 
These  improvements  increase  to  L=2S  days  where  improvements  of 
tne  order  of  20%  for  the  constant  drag  coefficients  and  30%  for 
the  the  linear  drag  coefficients  values  are  common. 

Figure  18(d)  shews  the  final  correlation  coefficients. 
Because  the  raw  correlation  coefficients  are  all  quite  high 
(>0.90),  the  broad  reatures  of  19(d)  are  not  significantly 
different  from  the  raw  values.  The  high  initial  correlations 
account  for  the  large  reductions  an  the  residual  variances  tnat 
can  re  -hieved  with  cnly  a  linear  regression.  Figure  19(d), 
nowevei ,  demonstrates  the  improvements  achieved  over  the  raw 
correlation  coe f f ic ien ts .  In  all  cases  the  improvements  are 
small  for  low  averaging  periods  and  rise  to  values  of  up  to  0.05 
at  1=28  days.  An  examination  of  Appendix  F  reveals  that  at  L=2 
days  there  are  no  values  less  than  0.99  .  Since  the  correlation 
ccefricient  is  rounded  by  1.0,  there  is  little  room  for 


i  rap  rove  men  t  at  the  lew  averaging  periods. 

For  averaging  periods  of  greater  than  1.0  days,  the  VA 
estimates  cf  the  long-term  variances  as  quantified  by  the  DV  and 
the  pci nt-fcy-pcint  variance  as  quantified  by  the  RV  can  be 
substantially  improved  through  application  of  the  £ 
corrections  given  in  Appendix  C.  In  regions  where  the 

unimproved  estimates  of  the  climatological  VA  fluxes  differed 
greatly  frem  the  1H  fluxes  (as  quantified  by  the  D Ms) ,  dramatic 
improvements  resulted  with  application  of  the  factor.  The 

major  strength  of  the  multivariate  regression  is  to  increase  the 
correlation  coefficients  particularly  at  longer  averaging 
periods.  If  a  similar  regression  had  been  performed 
irrespective  of  wind  speed  grouping  no  improvements  would  have 
resulted  because  the  correlation  coefficient  is  independent  of 
an  overall  scale  change. 
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CHAPTER  IV  HEAT  fLUXES 

This  chapter  describes  the  effects  of  averaging  winds, 
temperatures,  and  absolute  humidities  on  the  sensible  and  latent 
heat  fluxes.  lhe  sensible  heat  flux  is  a  measure  of  the  heat 
transferred  by  convection  into  or  out  of  the  sea  surface  while 
the  latent  heat  flux  is  a  measure  of  the  heat  lost  or  gained  by 
the  evaporation  cr  condensation.  The  bulk  aerodynamic  formulae 
for  their  calculation  are  given  in  Eguations  1.2(c)  and  (d)  . 

4. 1  Three-hourly  Heat  Fluxes 

Similar  to  the  3H  momentum  fluxes,  the  3H  heat  fluxes  had 
summer  minima  and  winter  maxima  average  values  of  which  station 
C  was  chosen  as  an  example. 

Plots  far  Station  C  of  the  yearly  and  monthly  latent  and 
sensible  3H  heat  fluxes  with  error  bars  indicating  t  1  sample 
standard  deviaticn  appear  in  Figures  20  (a)  and  (b) .  Similar 
graphs  with  ♦  1  standard  deviation  of  the  monthly  means  appear 
in  fiqures  20(c)  and  (d)  respectively.  In  this  section,  the 
sensible  heat.  results  will  be  quoted  in  parentheses  following 
the  latent  heat  results.  The  extreme  values  in  the  means  are 
much  greater  between  months  than  between  years.  Between  months 
the  maximum  values  are  8C.8  (34.7)  Watts/m2  in  December ( January) 
and  the  minimum  values  are  12.  4  (—  7. 4)  Watts/m2  in  July(July). 
The  minus  sign  indicates  a  net  transport  of  heat  from  the  air  to 
the  sea.  This  may  be  due  to  advection  of  warm  continental  air 
in  the  summer  to  the  mi d- A 1 1  ant ic.  Between  years  tbe 
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Figure  20.  due  year-to-year  and  monthly  variations  of  the  aean 
heat  fluxes  at  Station  c.  The  error  tars  are  t  1  sample 
standard  deviation  in  Figures  20(a)  and  (b)  and  ±  1  standard 
deviation  of  the  monthly  means  in  Figures  20(c)  and  (d). 


latent  (sensible)  heat  fluxes  have  maxiaa  in  1959  (1957)  of 
67.9(19.9)  Hatts/m2  and  rainiaa  in  1949  (  1949)  of  34.5(2.3) 
Hatts/a2. 

It  is  clear  froa  Figure  20  that  the  standard  deviations  in 
the  yearly  heat  values  are  fairly  constant.  The  aaxinua  yearly 
standard  deviations  for  the  latent  (sensible)  heat  fluxes  are 
90.  3  (69.8)  »atts/az  for  1957  (  1957)  while  the  niniuua  values  are 
62.  2  (34.4)  Hatts/a2  in  195o(195Q).  The  standard  deviations  of 
the  aonthly  leans  for  both  heat  fluxes  are  approximately  one- 
half  of  the  total  standard  deviations  based  on  all  the  three- 
hourly  flux  readings  indicating  that  about  one-guarter  of  the 
total  variance  nay  be  in  periods  of  greater  than  one  aonth. 

Figures  20  (t)  and  (d)  show  the  difference  in  standard 
deviation  between  +;he  suiter  and  winter  months.  The  Deceaber 
heat  flux  standard  deviation  is  91.1(b9.5)  Watts/a2  while  the 
July  values  are  38.  3(19. C)  Watts/m2.  The  standard  deviations  of 
the  monthly  means  vary  froa  34.0  (26.  8)  Watts/m2  in  Deceaber  to 
12.6(6.3)  Watts/s2  in  July  indicating  that  toe  annual  cycle  is 
present  in  the  Beans  and  the  variances. 

The  total  ships'  Beans  and  standard  deviations  of  the 
velocity  coaponents,  air-sea  teaperature,  air-sea  absolute 
humidity  difference,  latent  heat  and  sensible  fluxes  appear  in 
Table  X.  The  aean  heat  fluxes  are  all  of  the  saae  order  as  the 
standard  deviations. 

oceanographically,  the  data  Bay  be  divided  into  four  Bain 
tegiaes.  Stations  A,  E  ,  c,  and  8  are  in  Sub-Arctic  regions. 
They  are  characterized  by  sea  surface  teaperatures  of  about  1®  C 
and  colder  air  teaperatures.  As  one  proceeds  south  the  air 
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TABLE  X 

The  itean  and  standard  deviations  of  the  x  and  y  component  wind 
velocities,  the  air-sea  temperature  differences,  the  air-sea 
humidity  differences  and  the  latent  and  sensible  heat  fluxes. 
The  first  line  for  each  ship  is  the  mean  value  and  the  second 
line  is  standard  deviation. 
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temperature  increases,  decreasing  the  air-sea  temperature 
dif  ference. 

Station  D  lies  at  the  transition  region  between  the  Gulf 
Stream  and  North  Atlantic  Drift.  Here  the  water  is  relatively 
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warm  (about  5-10<>C)  accounting  for  large  air-sea  temperature 
differences.  The  advecticn  of  large  horizontal  sea  surface 
temperature  gradients  and  invasions  of  continental  air  account 
for  large  standard  deviations  in  the  heat  fluxes. 

Stations  E  and  N  lie  in  the  Kiddle  of  the  large  Atlantic 
and  Pacific  sub-tropic  gyres.  They  are  characterized  by 
Moderate  air-sea  temperature  differences  and  by  dry  prevailing 
winds  thus  accounting  tor  noderate  sensible  heats  and  large 
latent  heat  flux  values.  Station  N  is  the  Trade  Hind  region, 
there  steady  summer  winds  account  for  the  small  heat  flux 
standard  deviations  in  comparison  to  Station  E. 

Stations  I,  J,  and  K  lie  on  the  eastern  boundary  of  the 
Atlantic  Ccean.  Proceeding  south,  the  air  temperature  becomes 
warmer  more  quickly  than  the  corresponding  increase  in  sea 
surface  temperature  accounting  for  lower  sensible  heat  fluxes. 
The  three  stations  show  similar  air-sea  humidity  differences  and 
moderate  latent  heat-  flux  values. 

4.2  Uncorrected  Test  He suits 

The  difference  means,  difference  variances,  residual 
variances  and  correlation  coefficients  as  defined  in  Equations 
2.16-2.19  were  calculated  with  no  correction  factors  applied. 
The  results  appear  in  Figu  res  21(a)  through  (h)  while  tne  raw 
values  are  in  Appendix  E.  Table  XI  contains  the  3H  heat  flux 
variance  at  each  averaging  period  so  the  difference  variances 
and  residual  variances  can  be  converted  to  absolute  quantities. 

As  in  the  wind  stress  case,  the  heat  flux  difference  means 
(DHs)  appear  to  te  dependent  upon  the  average  3H  heat  flux  value 
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TABLE  XI 

The  absolute  latent  and  sensible  heat  flux  variances  as  a 
function  of  averaging  period.  The  columns  refer  to  the 
averaging  period  in  davs  and  the  variances  are  in 
(Watts/m*)  *x10 *. 


|  VARIANCES  -  SENSIBLE  HEAT  FL0I  J 

|  SHIP  |  0.25  |  0.50  |  1.0  |  2.0  |  4.0  j  7.0  |  14.0  J  28.  0  | 

|  A  |  3.511|  3.  303  J  2.  9821  2.  508)  2.035)  1.669}  1.331J  0.994J 

|  E  |  6.154|  5.  900|  5.461|  4.  870)  4 .  1 4  2 1  3.6841  3.077|  2.5301 

|  C  |  2.496|  2.314)  2.C201  1.626|  1.2631  1.045}  0.797}  0.566J 

|  D  |  6.947J  6.512|  5.7251  4.6461  3.6791  3.084|  2.4791  1.892J 

!  E  |  1.  1  331  1.  05  7  }  0.921|  C.  759J  0.573!  0.451J  0.  3261  0.245} 

}  I  |  2.174}  2.219|  1.998|  1 .  607 1  1.38J|  1.  1331  0.8671  0.670| 

|  J  |  1.423|  1.316)  1.160|  0.968)  0.747}  0.6021  0.412)  0.3251 

|  K  |  1.0011  0.941)  0.860|  0.764J  0.672|  0.589)  0.510J  0.442| 

1  M  |  2.398|  2.266)  2.058)  1.775)  1.453)  1.200)  0.900J  0.7C6J 

|  N  |  0.  341)  0.  302|  0.  2751  0.  23b)  0.  1 92 1  0.1b2|  0.  124)  0.081| 

|  F  )  0.8561  0.  806)  0.7C8)  0.  577|  0.4481  0.3451  0.240)  0.  181J 

I  I  I  I  I  I  1  1  I  1 

} - L - JL - L - X _ I - 1 _ I _ i _ 1 

)  VAFIANCES  -  LATENT  H  BAT  FLUX  1 

1 - + - j. - + - -} - f - ) - + - j - 1 

JSHIF)  0.25  |  0.50  |  1.0  |  2.0  |  4.0  |  7.0  1  14.0  1  28.0  ) 

l - i - + - + - H - + - d - 1 - i- - 1 

|  A  |  4.  078)  3.  785)  3.  343)  2.  751  J  2 .  179)  1 . 8 12 J  1.426J  1.064) 

|  E  |  4.017|  3.795|  3.443)  2.989)  2.4841  2.161|  1.812)  1.518| 

|  C  J  5.326)  4.951)  4.3301  3.495)  2.663)  2.1541  1.639}  1.233} 

|  D  1  21.51}  20.  15)  17.70)  14.71  |  11.67|  9.680|  7.830}  6.171} 

)  E  |  12.79)  12.01)  10.74)  9.036)  7.008)  5.b93)  4.237}  3.295} 

|  I  |  4.600)  4.  205)  3.  655}  2.  933  }  2.318}  1.879)  1.373)  1.029| 

1  J  |  5.0  10)  4.  582)  3 . 96b  j  3.  1  12)  2.242)  1.732)  1 .  1 32  J  0.836) 

|  K  |  6.  155)  5.  753)  5.  1651  4.  351)  J.500|  2.8471  2.0911  1.668J 

I  B  I  3.374J  3.122!  2.747)  2.297|  1.838)  1.502}  1.133}  0.884} 

J  N  |  5.084)  4.701)  4.  2  13 J  3.  542|  2.7321  2.  1571  1.481}  0.870| 

|  F  |  3.11b|  2.653)  2. 509 |  2.034)  1.580}  1.2531  0.957}  0.750} 

I  I  I  I  I  I  1  I  I  1 


and  with  minor  exceptions,  increase  approximately  exponentially 
with  averaging  period.  A  comparison  jf  Table  X  and  Appendix  E 
indicates  that  the  sizes  of  the  difference  means  are  well 
correlated  with  the  average  3H  heat  flux.  At  L=28.0  days,  the 
sensible  heat  fluxes  fall  into  two  distinct  groups  --  Stations 
A,  E,  D,  1,  and  r  with  average  DM s  greater  than  18  Watts/m*  and 
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Stations  C,D,  E,  J,  N,  and  P  having  QMs  less  than  10  iatts/a*. 
The  latent  heat  flaxes  have  three  distinct  groups.  Stations  D 
and  E  have  uncorrected  DMs  greater  than  60  Watts/*2.  At  Lx28.0 
days.  Stations  A,  B,  C,  I,  J,  K,  M  and  N  have  DMs  in  the  range 
of  20  to  60  Watts/o2;  Station  P  has  distinctly  saaller  DMs  at 
less  than  20.0  Vatts/a2.  The  saae  groupings  occur  in  the  3H 
aean  flux  values  shown  in  TaDle  X  indicating  that  the  larger 
cli vatological  errors  occur  in  regions  with  the  higher  initial 
tean  flux. 

The  difference  variances  (DVs)  for  all  ships  excluding 
Station  N  differ  only  slightly.  Both  heat  fluxes  show  increases 
fron  near  0  at  1=0.26  days  to  mean  Baxiaun  values  of  about  0.55 
to  0.65  at  L=28. 0  days.  Station  N  has  aarkedly  lower  DVs  from 
1=1.0  to  L-14.0  days  in  the  sensible  heat  flux.  For  the  latent 
heat  flux,  stations  N*s  values  are  nearly  0.0  Watts/n2  up  to 
1=4. C  days  and  then  go  negative  --  in  contrast  to  all  the  other 
ships. 

All  ships  have  siailar  residual  variances  (HVs)  up  to  L=4 
days.  The  sensible  heat  fluxes  then  split  into  two  distinct 
groups.  The  first  includes  Stations  A,  E,  1,  A,  and  N  having  BV 
values  of  qreater  than  40. 0T  at  L=28.0  days.  The  others  have 
F  Vs  less  than  30. OX.  The  latent  heat  flux  BVs  are  generally 
larqer  (except  at  Station  P)  than  for  the  sensible  heat  fluxes. 
The  averaqe  value  at  L=28.C  days  is  about  35*.  Up  to  L=7.0 
days,  the  Station  N  values  are  indistinguishable  (although 
biased  hiqh)  frc»  the  other  values.  At  L=14.0  days,  however, 
the  Station  N  values  are  such  greater  than  the  other  ships. 
Indeed,  at  L=  28.0  days  the  Station  N  RV  is  greater  than  100.0* 
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which  is  net  shewn  in  Figure  21. 

The  correlation  coefficients  for  the  heat  fluxes  are  lower 
than  for  the  momentum  fluxes.  Those  for  the  latent  heat  flux 
are  lower  than  for  the  sensible  heat  flux.  At  L=28.0  days 
Station  has  the  lowest  correlation  coefficient  of  0.438  (not 
shown)  for  the  latent  heat  flux  and  0.736  for  the  sensible  heat 
flux.  After  L=1h.O  days,  the  sensible  heat  flux  correlation 
coefficients  divide  into  the  sane  two  groups  as  for  the  fiV 
values  --  Stations  A,  E,  I,  H ,  and  N  with  values  less  than  0.835 
(at  1-28.0  days)  and  the  rest  having  correlation  coefficients 
greater  than  C.900. 

4.3  H eajt  Regression 

Figure  22  shows  an  example  of  the  3H  heat  flux  plotted 
against  the  VA  heat  flux  for  Station  A,  1=28  days.  Siailar  to 
the  stress  case,  shown  in  Figure  6,  the  scatter  increases  for 
ex*reme  positive  and  negative  values.  As  the  VA  flux  approaches 
0.0  Watts/m2,  the  scatter  becoies  wore  constricted.  Mote, 
however,  that  there  is  a  bias  towards  positive  heat  fluxes. 

Twc  important  differences  exist  between  the  heat  flux 
regression  and  the  momentum  flux  regression.  First,  the 
analysis  of  the  raw  3H  values  revealed  that  the  latent  and 
sensible  heat  fluxes  are  statistically  guite  distinct 
quantities.  Consequently  the  regression  was  perforaed 
separately  cn  each  flux.  Thus,  no  attempt  was  aade  to  aasuae 
that  the  fluxes  siqht  be  similar  and  thus  produce  a  single 
t r  insfor mat  ion  applicable  to  both.  Second,  as  shown  in  Table  X, 
the  means  art  of  the  save  aagnitude  as  the  standard  deviations. 
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Figure  22.  The  sensible  and  latent  2ii  versus  VA  heat  fluxes  at 
Station  C ,  L=23.0  days. 


consequently  the  following  quantity: 
1 3  X**. 
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WdS  minimized  as  outlined  in  Appendix  A.  1  where  X  and  X'  are  the 
jn  and  VA  variates  f c:  either  tn^  latent  or  sensible  neat 
rluxts.  This  gives: 
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ijuation  of  the  form 


-  b 


+  ■  v  x 


wnere  x  and  b  are  vectors  with  components  corresponding  to  13 
Eeaufort  categories  and  _A_  is  a  13x13  matrix.  The  solution  was 
obtained  by  Gaussian  elimination  followed  by  iterative 
improvement  until  a  relative  error  of  less  than  .001  was  found. 

Appendix  A.  1  shows  that  the  DVs  were  exactly  equal  to  the 
AVs  ar.d  the  systematic  positive  DV  bias  again  exists.  This  fact 
provided  a  convenient  check  to  ensure  that  "-he  computer 
programmes  were  running  correctly.  The  systematic  bias  was 
eliminated  by  : 


—  •  w  k  -e 
'•IV 


wnero  v'\i  is  the  3i!  standard  deviation  ,  'd’V 1  is  the  VA  standard 
deviation  with  ice  applied  and  ^  represents  either  the  latent 
or  sensible  heat.  An  attempt  was  made  to  correct  on  a 


tieauiort  category  rasis  (ie.  %*>.<.  ).  However, 

<5Vyi_ 

cases  the  requited  iterations  failed  to  converge. 


in  many 


Anpenlix  A.  3  shows  that  wh*^n  readjusting  from  the  minimum 
RV  value,  the  correction  given  in  Equation  4.4  induces  a  further 
error  i n  R  V  given  by: 
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error  is  less  than  ,.C,  . 

Again,  the  primes  will  be  dropped  and  it  will  be  understood 
that  includes  the  DV  correction. 

4.4  Wind  Dependent  Correction  Factors 

The  £««for  the  heat  fluxes  were  calculated  and  appear  in 
Appendix  F  -  Part  I.  Figure  23  shows  the  values  as 

functions  of  wind  speed  and  averaging  period. 

In  general,  the  inverse  heat  corrections  follow  a 

similar  pattern  to  the  stress  values.  For  low  wind  speeds 
approaches  0.0  (  lh.eapproaches  co  )  .  As  the  wind  speed 

increases,  approaches  1.0  for  all  averaging  periods.  In 

general  Appendix  F  indicates  that  the  corrections  required  for 
tne  sensible  heat  fluxes  are  larger  than  those  required  for  the 
latent  heat  fluxes.  Furthermore,  a  comparison  with  Appendix  C 
indicates  that  the  stress  corrections  are  larger  than  both  the 
heat  corrections. 

An  estimate  cf  the  error  in  £ was  obtained  by  employing 
the  Student  't1  test  as  outlined  in  Equation  3.10.  The  results 
appear  in  Appendix  F  -  Part  II.  ihese  results  should  be  viewed 
somewhat  skeptically  as  no  attempt  was  made  to  assure  that  the 
scatter  about  the  regression  was  Gaussian.  Similar  patterns  to 
the  stress  errors,  however,  are  exhibited  in  the  heat  flux 
errors.  For  the  lowest  wind  speed  categories,  the  95% 

confidence  zone  is  large  compared  to  the  £ ^  values.  By 
Beaufort  interval  3  (1.6  -  3.4  m/sec2)  errors  in  the 

calculations  are  reduced  to  less  than  10%  cf  the  calculated 
value  for  all  averaging  periods.  In  general,  the  errors 
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increase  as  the  averaging  period  increases. 

The  £yi t  heat  corrections  for  all  ships  are  fairly  tigntly 
grouped.  Station  N  shows  significantly  larger  1.0/  (smaller 

)  values  at  averaging  periods  greater  than  0.5  days. 
Station  K  also  shows  significantly  larger  1.0/^^  values  for  the 
sensible  heat  flux  —  lying  between  the  Station  N  values  and  the 
main  ship  grouping.  An  inspection  cf  the  sensible  heat  raw 
statistics  in  Table  X  does  not  reveal  anything  either  to 
distinguish  Station  K  from  the  other  ships  cr  to  show  any 
consistencies  with  Station  N  to  account  for  the  difference.  In 
the  latent  heat  fluxes,  the  Station  K  values  are  biased  towards 
high  values  in  the  plots  tut  can  be  considered  within 

tne  general  ships'  geographic  grouping  and  are  far  below  those 
of  Station  N. 

4.5  Beaufort  Grouted  Test  Results 

The  three  test  quantities  were  calculated  for  the  iT*.* 
values  given  in  Appendix  F.  The  results  appear  in  Appendix  G 
and  are  shewn  in  Figure  24.  The  arithmetic  differences  between 
Figures  24  and  21  are  shown  in  Figure  25.  The  difference 
variances  are  not  considered  in  this  section  because  the 
correction  for  the  systematic  bias  forced  all  values  to  exactly 
0.0. 

Tne  difference  means  have  improved  markedly  after  the 
application  of  A  comparison  of  Figures  24  and  21  indicates 

that  noth  the  latent  and  sensitle  heat  flux  DMs  have  been 
improved  ty  a  factor  cf  5.0.  Climatological  errors  are  now  less 
tnan  3.5  Hatts/m*  for  the  sensible  heat  flux  and  17.7  Watts/m* 
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Figure  24.  lest  values  with  individual  ships* 
the  heat  fluxes. 
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Figure  25.  Test  value  improvements  over  Figure  21,  t 
uncorrected  case. 


for  tne  latent  heat  flux.  Ignoring  Stations  E  and  N  for  the 
latent  heat  flux,  the  DM  errors  are  then  all  less  than  10.0 
Watts/m2.  Ihe  improvements  shown  in  Figure  25  and  are  virtually 
a  replication  of  Figure  21  indicating  that  the  greater  the 
initial  DM  error,  the  greater  the  improvements  that  can  be 
obtained. 

The  residual  variance  errors  are  somewhat  less  encouraging 
for  longer  averaging  periods.  Station  N  in  particular,  has 
extremely  high  (up  tc  0.23  for  the  sensible  heat  and  0.83  for 
the  latent  heat  flux)  RVs  at  L=14.0  and  L=28.  0  days.  If  one 
considers  10%  to  be  an  acceptable  noise  level  then  an 
appreciable  number  of  ships  lie  above  this  level  at  L=14.0  days 
for  the  sensible  heat  flux  and  at  L=4.0  days  for  the  latent  heat 
flux.  The  application  of  £1^  has  markedly  improved  the  final 
residual  variances  over  the  initial  raw  values.  The  RVs  at  long 
averaging  periods  are  now  about  one-half  the  uncorrected  values. 
Figure  25  indicates  that  improvements  of  0.25  for  the  sensible 
neat  rlux  and  0.15  for  the  latent  heat  flux  at  L=28.0  days  are 
typical. 

The  application  of  £  (*.4  has  improved  the  correlation 
coefficients.  The  sensible  heat  flux  correlations,  with  the 
exception  of  Station  N,  L=28.0  days,  are  all  above  0.90.  The 
latent  heat  flux  coi  .relations  are  above  0.90  up  to  L=4.0  days 
after  which  seme  stations  fall  off,  for  example,  to  values  as 
low  as  of  0.583  at  Station  N  and  0.750  at  Station  E  (both  not 
shewn)  for  L=28.Q  days.  According  to  Equation  4.8,  this  partly 
explains  the  large  RVs  at  Stations  N  and  E.  if  the  DVs  were  not 
corrected,  the  RVs  would  be  reduced  by  17.4%  and  6.25% 


respectively.  The  improvements  shewn  in  Figure  25  have  similar 
ranges  for  both  the  latent  and  sensitle  heat  fluxes.  There  is  a 
tendency  for  the  stations  with  the  lower  initial  correlation 
coefficients  to  shew  the  most  improvement.  For  example,  in 
Figure  21,  the  sensible  heat  flux  raw  correlations  at  Stations 
A,  E,  I,  H  and  B  are  distinctly  grouped  lower  than  the  others. 
Figure  25  shows  these  same  stations  grouped  with  the  most 
impro vement. 
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CH&PTEB  V  EMPIRICAL  FORMULA  AND  TEMPORAL  VARIATIONS 


5. 1  Introduction 

In  the  previous  two  chapters  it  has  been  established  that 
systematic  reductions  cccur  in  the  3H  momentum  and  heat  fluxes 
when  estimated  through  VA  parameters.  In  Chapter  III  it  was 
shown  that  the  ^/?^l  optimally  regressed  stress  reductions  were 
qualitatively  similar  tc  the  set  of  averaged  corrections  R  (L) . 
With  the  exception  of  Station  K,  sensible  heat  flux,  it  was 
noted  that  the  1/£*.t  optimal  corrections  for  the  heat  fluxes 
were  similar  tc  the  corresponding  l/t*.*.  corrections  found  in  the 
stress.  R j  (L)  of  Equation  3.14  can  te  rewritten  as: 


/r.U'i  <  1  *  AieL 


5.  1 


using  notation  previously  defined.  As  stated  previously  the 
variance  information  is  lost  to  the  averaging  process,  however, 
the  VA  wind  speeds,  (u;2*  v2)  i/2 ,  are  required  for  the  estimates 
of  all  the  VA  fluxes.  Given  these  arguments,  an  expected  wind 
speed  dependence  of  £>«».  for  all  fluxes  is: 


t 


tec 


1  +  o<  (  V 


5.  2 


where  and  ^  are  constants  to  te  determined. 

The  temporal  aspects  of  the  reductions  have  also  not  been 


investigated. 

The  se 

may  manifest  themselves  in 

two  ways. 

First,  there  may 

be 

consiste  nt 

variations  between 

averaging 

periods  which 

will 

te  called 

functional  temporal 

variations. 
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Establishing  consistencies  in  these  variations  would  aid  greatly 
in  interpolating  the  reductions  to  averaging  periods  not 
specifically  investigated.  Second,  the  ^lu.  transformations  may 
be  influenced  by  time  scales  greater  than  the  averaging  period, 
for  example,  the  corrections  required  in  December  may  be 
markedly  different  from  those  required  in  July.  These  will  be 
called  intrinsic  temporal  variations. 

Finally,  it  was  shown  in  Chapters  III  and  IV  that  the 
optimal  corrections  (with  the  exception  of  Station  N  for  all 

fluxes,  and  Station  K  for  the  sensible  heat  flux)  dc  not  differ 
markedly  from  one  ship  location  to  the  next.  If  a  formula  which 
is  applicable  at  most  ship  locations  can  be  found,  then  the  ease 
of  applying  the  tranfcrmations  would  be  greatly  enhanced. 

5.2  Empirical  Formula 

To  verify  the  arguments  of  the  previous  section,  the 
values  of  Appendix  C  for  the  stresses  and  of  Appendix  F  for  the 
heat  fluxes  for  each  ship  were  weighted  and  averaged  (excluding 
Station  N)  according  to  the  method  outlined  in  Section  3.3  to 
arrive  at  a  geographically  averaged  estimate  3m.  of  the 
corrections.  The  results  appear  at  the  ends  of  Appendices  C  and 
F  -  Farts  I. 

(  fc*,,  "1)  was  initially  plotted  versus  (u**v*)  &  on  a  log 
scale  and  the  results  for  all  fluxes  appear  in  Figure  26.  The 
horizontal  line  indicates  where  a  correction  of  10X  or  greater 
is  required.  In  all  cases  the  curves  are  quite  linear  in  the 
range  0.5  to  20  m/sec  except  in  the  linear  drag  coefficient  case 
where  a  discontinuity  in  the  slope  occurs  at  10  m/sec.  This  is 
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Figure  26.  The  - 1  values  as  a  function  of  wind  speed.  All 
points  above  the  horizontal  line  reguire  greater  than  10% 
correction.  The  inset  numbers  indicate  averaging  period  in 
days. 


the  wind  speed  at  which  the  drag  coefficient  changes  form  (see 
Equation  2.5).  For  wind  speeds  greater  than  20  m/sec,  the 
points  neccme  scattered.  This  may  be  due  to  there  being 


substantially  fewer  averagings  in  this  regime.  The  slopes  for 

a lx  curves  appear  to  be  similar  for  all  averaging  periods.  Note 

that  typical  wind  speeds  of  1-10  m/sec  at  L=28.0  days  yield 

corrections  of  about  1  to  8  which  are  of  the  correct  order  to 

improve  the  transport  calculations  by  Aagaard  shown  in  Figure 

1(b).  The  individual  ships'  plots  of  (?r«.-1)  versus  ( u,2 *■  v-2 )  »/2 

0  * 

were  examined  and  were  all  qualitatively  similar  to  the 
geographically  averaged  case. 

Next,  (£>{-1)  was  plotted  versus  averaging  period  (L)  on  a 
log  scale  tc  determine  any  functional  temporal  variations  which 
appear  in  Figure  27.  Again  all  fluxes  show  similar  patterns. 
All  points  above  the  horizontal  line  require  at  least  10% 
correction.  Beaufort  interval  1  (0.2  m/sec)  is  more  scattered 
than  the  others.  This  extends  to  Beaufort  interval  2  (1.0 

m/sec)  on  the  individual  ships  (£«u-1)  versus  L  plots.  This  is 
not  surprisinq  since  the  extremities  in  wind  speed  contain  the 
fewesx  number  of  points  and  the  lowest  wind  speeds  have  the 

largest  inherent  errors  of  Appendices  C  and  F  -  Parts  II.  At 

L=2.0  days,  the  ^  represents  Beaufort  category  12  (26.0 
m/sec)  and  appears  distinct  from  the  pattern  of  the  rest  of  the 
points.  This  point  was  determined  from  only  four  out  of 
approximately  one  million  total  possible  averagings  and  cannot 
be  considered  statistically  valid. 

All  ether  curves  appear  to  increase  systematically.  From 
L=0.25  to  L=2.0  days  (which  I  call  Region  I)  the  slopes  are 

steeper  than  frem  L=4.0  tc  L=28.0  days  (which  I  call  Region  II). 
The  vertical  line  at  L=3.0  days  demarcates  the  two  regions. 
Within  each  region  the  slopes  do  not  appear  to  change 
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Figure  27.  Ihe  £*e-1  values  as  a  function  of  averaging  period. 
The  vertical  line  demarcates  the  division  between  Regions  I  and 
II  and  all  points  above  the  horizontal  line  require  greater  than 
10%  correction.  Ihe  inset  numbers  indicate  wind  speed  in  m/sec. 


dramatically  frcm  one  Eeaufort  category  to  the  next. 

Figure  26  indicates  a  power  law  behaviour  of  (^w-c-1)  with 
wind  speed  which  appears  to  be  independent  of  averaging  period 
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and  Figure  27  indicates  a  power  law  behaviour  of  (£**"1)  with 
averaging  period  which  appears  to  be  independent  of  wind  speed. 
It  is  assumed  a  priori  that  the  correction  7^  is  0f  the  form: 

,  v  &t  *■ 

1-  L 

wnere  (u f+vf)ll2  is  the  vector  averaged  wind  speed  and  L  is  the 
averaging  period  in  days. 

The  assumed  form  of  the  X  (3H)  variate  and  X'  (VA)  residual 

is  : 

g\  &x  =  Cx-  nx'f 

5.4 

where  the  cverbar  indicates  averaging  over  all  Beaufort 
categories  and  averaging  periods  which  can  be  rewritten  as: 

<r?s\=  l  u-r) >~)*T 

which  immediately  leads  tc  : 

<ri  C«-  Ux-vW^^VxT  5.5 

Note  tnat  we  new  regress  the  assumed  form  directly  to  the  data. 
This  is  a  ncn-linear  regression  of  the  form: 

<y'  s\  ■  (Y-  ?(x,  ;  L)  X'  )*■  5. 6 

where  Y  =  X  -  X',  X.'  =  o'  ,  A  ,  ^  for  i=1,  2,  3  and 
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tiki  (U-J2*  Vj2)  »/*  ,  L)  =  )PI*  LY  .  If  we  have  an  initial  guess 

f  (  Xt,  ),  Eguaticn  5.6  can  be  expanded  in  a  Taylor  series  to 
first  order  as  fellows: 


The 


quantity  AA  ^  can  then  be  minimized  giving: 
0=  - 


5.8 


this  is  a  matrix  equation  of  the  form: 


r  -  A* 

where: 

rr 

Am-  "Si  SAA 

DV.|Xc0 

Kl  -  A  k 


5.9 


1 
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by  evaluating  explicitly  the  derivatives  in  Equation  5.8,  the 
can  then  be  iterated  until  AV  is  snail  compared  to 

.  Iteration  was  halted  when  I /Xio\  ^  o  ° 0  1  .  Note 

that  XX'  and  X'2  are  the  cross  products  and  squares  of  the 

individual  variates  and  not  the  squares  of  the  average  within 
each  beaufort  category.  Furthermore,  the  difference  means  have 

oeen  removed  from  the  regression  residuals.  It  will  be 

demonstrated  that  this  has  little  influence  on  the  final  result. 

Since  Figure  27  indicated  a  definite  break  in  the  slope  at 
L=3  days,  the  regression  was  performed  piecewise  with  a  separate 
set  of  constants  for  tie  two  regions.  A  regression  of  this  type 
greatly  reduces  the  required  number  of  constants  to  be  stored; 
consequently  a  separate  regression  was  performed  on  the  x  and  y 
stress  components.  Thus  the  number  of  parameters  were  reduced 
from  104  per  ship  to  12  for  the  stresses  and  6  for  each  of  the 
heat  fluxes. 

To  be  rigorous,  the  regression  should  be  performed  directly 
against  the  VA  and  3H  estimates  from  each  individual  averaging. 
Tnis  would  entail  calculations  involving  as  many  as  50,000 
points  per  snip  and  flux  repeated  over  several  iterations.  The 
computer  costs  would  be  excessive  in  storage  and  processing 
time.  The  Eeauf ort-a veraging  period  groupings  of  the  sums,  sums 
of  squares  and  sums  cf  cross  products  of  the  3H  and  VA  variates 
were  stored  during  the  calculations.  By  appropriately 

weighting  these  stored  statistics,  an  approximation  for  the  , 
^  and  X  could  be  achieved  requiring  effectively  only  13x8*104 
feints  per  ship  and  flux. 

A  weighting  matrix  was  created  in  order  not  to  bias  the 
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parameter  estimates.  First,  each  Beaufort  category  contained 
unequal  numbers  of  points.  The  first  order  weight  then  was 
directly  proportional  to  the  number  of  3H/VA  averagings  within 
each  Beaufort  category.  By  the  nature  of  the  averaging  of  the 
3H  and  VA  variates,  the  longer  averaging  periods  had  fewer 
estimates  of  the  reduction  which  would  tend  bias  the  temporal 
power  (  t  )  in  favour  of  the  shorter  averaging  periods. 
However,  we  wish  an  unbiased  estimate  over  all  averaging 
periods.  Consequently,  the  number  of  estimates  within  each 
Beaufort  category  and  averaging  period  was  expanded  to  the  total 
number  of  raw  three-hourly  samples  that  went  into  the 
calculations  of  the  statistics  of  each  category.  Thus,  if 
Beaufort  interval  8  at  L=1.0  days  contained  4  3H/VA  estimates, 
since  there  were  8  readings  per  day,  the  4  readings  were 
expanded  to  32. 

This  type  of  regression  minimizes  the  total  residuals 
within  each  region.  That  is,  each  averaging  period  is  no  longer 
considered  as  a  single  block  of  data  but  the  total  residuals 
from  1=0.25  to  L=2. 0  days  (Region  I)  and  from  L=  4.0  to  L=28.0 
(Region  II)  are  now  minimized.  It  will  be  shown  that  the 
regression  performed  in  this  manner  did  not  significantly 
increase  the  residuals  when  compared  with  a  direct  regression. 

The  first  guess  {  \Vo  )  was  made  by  performing  a  linear 
regression  in  log  space.  At  several  locations,  the  iteration 
failed  to  converge  with  this  value.  Fortunately,  it  was 
established  that  most  ships  converged  to  values  of  {Q  =-1.0  and 
^  =1.0  in  Region  I  and  V  =0.25  in  Region  II.  Consequently 
these  values  as  well  as  an  arbitrary  oi  =2.0  were  used  as 
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initial  guesses.  No*  only  did  all  ships  converge,  bat  in  cases 
where  the  iteration  had  previously  succeeded,  the  convergence 
vent  to  the  same  values  in  fever  recursions.  This  indicates 
that  a  linear  regression  in  log  space  nay  not  produce  suitable 
estimates  of  the  coefficients. 


5.3  Ship  Paranetei  Estimates 

In  light  cf  the  draaatic  improvements  in  test  results 
achieved  by  the  calculations,  they  provided  an  acceptable 

means  of  comparison  to  determine  the  effectiveness  of  the  ?7 
technigue.  If  t  he  V  results  compare  favourably  with  the 
results,  then  the  empirical  formula  would  be  an  effective 
replacement.  In  this  section,  the  residual  variances  mere 
calculated  by  applying  the  regression  coefficients  against  the 
stored  Beaufort  number-averaging  period  statistics  outlined 
above  and  the  empirical  formula  was  not  allowed  to  interpolate 
wind  speed. 

Pigure  28  shows  the  "Y  FV  estimate  minus  the  BV 

estimate  in  percent.  For  the  stresses,  nearly  all  results  are 
within  1.59  of  the  values  for  all  averaging  periods.  The 

only  exception  is  Station  D,  L=7.0  days,  linear  drag 
coefficient,  x  component  where  the  difference  mas  2. 3f  (not 
shown).  Station  N  was  also  within  1.5%  of  the  £*<  residuals. 
Many  stations  showed  an  improvement  (ie.  negative  residual 
variances  in  Figure  28)  in  the  stresses.  This  occurred  because 
here  a  separate  regression  was  performed  on  the  components. 

The  heat  flux  RVs  (Figures  28(e)  and  (f))  were  generally 
within  4.0%  of  the  RVs  for  all  averaging  periods.  The 
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exceptions  were  Station  N  (not  shown)  where  the  results  in 
Region  II  were  generally  poor  and  Station  J,  L- 26,0  days,  latent 
heat  flux  where  the  regression  induced  a  further  error  of  about 
10. OX  (also  not  shown).  Note,  however,  up  to  L=14.0  days  the 
differences  were  all  within  1. 5%  of  the  values. 

As  a  further  indication  of  the  accuracy  of  the  technique, 
ali  fluxes  (momentum  and  heat)  showed  a  positive  bias  in  the 
difference  variances  (not  shown).  Although  they  were  not 
identical  to  the  residual  variances,  they  were  close  enough  to 
indicate  that  a  DV/RV  statistical  biasing  was  occurring, 
consequently,  oi.  cf  Equation  5.3  was  adjusted  by  a 
multiplicative  factor  to  reduce  the  difference  variances. 
Again  considering  the  3H  variate  to  be  X  and  the  VA  variate  to 
be  X',  the  condition  for  0  difference  variance  for  each 
averaging  period  is  : 


where  the  summation  is  over  13  Beaufort  intervals 


.  ¥ 


the  minimum  approximation  outlined  in  Equation  5.6,  (4  is  a 

multiplicative  correction  factor  forcing  the  difference  variance 
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tc  0  for  each  averaging  period,  and  X'nt=2Xiwv  •  This 
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where  ^  ^  T are  the  raw  3H  and  VA  variances  with  no 

correcticn  applied.  This  is  then  a  quadratic  in  Af  •  The 
correct  root  was  chosen  by  comparison  with  the  uncorrected 
difference  eean  value.  The  corrections  for  each  averaging 
pericd  within  each  region  were  then  averaged  to  detereine  a  eean 
regicnal  correction. 

After  this  correction  was  applied,  the  constants  *  •  ft  • 

and  ^  of  Equation  5.3  were  calculated  and  the  results  appear  in 
Appendix  H.  The  ^  tera  for  the  linear  drag  coefficient 
0.9),  constant  drag  coef f icient  (  ~  -1 . 3) ,  and  heat  fluxes  (-1.2 
to  -1.4)  tend  to  decrease  in  this  order,  for  aost  ships.  The 
teapcral  power  (  ^5  coefficient)  is  regular  for  all  ships 

excluding  Station  N.  In  Region  I,  it  varies  from  a  low  of  0.860 
to  a  high  of  1.04  indicating  that  the  corrections  required  in 
this  region  vary  approximately  linearly  with  averaging  period. 
In  Begion  II  the  ^  coefficients  vary  from  0.16b  to  0.359  with 
aost  lying  between  0.20  and  0.30.  Furtheraore  the  temporal 
power  is  independent  of  type  of  flux.  Thus  the  discontinuity  in 
the  tine  dependence  is  quite  distinct  between  two  and  four  days 
which  appears  to  correspond  to  a  discontinuity  observed  in  the 
slope  of  the  stress  coaponents  spectra  noted  by  Hillebrand 
(1978)  . 

The  four  test  functions  were  then  calculated.  Calculation 
of  the  values  required  storing  the  sums,  suas  of  squares, 

and  suas  of  the  cross  products  of  the  38  and  VA  variates  within 
each  Beaufort  category.  These  stored  values  served  as  the  raw 
statistics  frca  which  the  four  test  functions  were  calculated. 
The  empirical  formula  allowed  interpolation  of  the  VA  wxnd  speed 
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correction  within  each  Beaufort  category.  (then  the  empirical 
formula  was  applied  to  each  ?a  estimate  individually,  certain 
ships  showed  a  dramatic  decrease  in  residual  variance.  For 
example.  Station  N  latent  heat  flux,  had  a  R?  of  82. OX  for 

1=28  days.  8hen  the  eapirical  foraula  interpolated  the 
velocities,  the  RV  was  only  58. 9%.  In  the  sensible  heat  flux, 
the  change  was  froa  23. 3X  to  12. 6X.  At  L=28.Q  days,  78X  of  the 
data  is  qrcuped  in  Beaufort  categories  3  and  4.  Consequently 
the  interpolation  of  the  correction  within  the  two  groups  was 
important.  At  locations  where  the  data  were  sore  evenly 
distributed  auong  Beaufort  categories  (eg.  Station  C) ,  the 
improvement  by  direct  application  of  was  less  than  1.0X. 

This  immediately  suggests  that  the  empirical  formula 
regression  coefficients  should  also  be  calculated  directly 
against  the  data  rather  than  against  the  Beaufort  grouped 
statistics.  Since  the  regression  involved  several  iterations  of 
large  amounts  of  data  the  only  feasible  method  was  to  regress 
the  data  in  Reqion  II  only.  This  was  done  for  several  ships 
including  Station  N,  latent  heat  flux  and  in  no  case  were  the 
residuals  improved  more  than  1.0X.  In  some  cases  the  difference 
variances  increased  to  values  of  about  30. OX.  Since  the 
regression  against  the  data  required  several  iterations  over 
large  amounts  cf  information  proving  extremely  costly  in 
computer  tine  and  was  of  dubious  value  ,  this  avenue  was  not 
pursued. 

The  four  test  functions  were  then  calculated  by  applying 
H  f o  each  VA  flux  and  stress  estimate.  The  complete  results 
appear  in  Appendix  I  while  the  x  component  linear  drag 
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coefficient  and  the  heat  fluxes  appear  in  Figures  29  to  32. 

The  difference  leans  appear  in  Figure  29.  Except  for 
Station  A,  all  DMs  are  below  0.08  dPa  for  the  linear  drag 
coefficient.  A  comparison  of  Appendix  I  with  Appendix  D  ( 


results) 

indicates 

than 

the  values 

found  by 

the 

empirical 

formula. 

although 

net 

identical  to 

the  £ « 

DMs, 

noth  show 

similar  ranges  of  values.  The  heat  flux  DMs  are  generally  1.0- 
2.0  Watts/m2  larger  with  the  empirical  formula  than  with  the 
directly  applied  technique.  The  sensible  heat  flux  DMs, 
with  the  exception  of  Station  I,  L=28.0  days,  are  all  less  than 
2.5  Watts/m2  while  the  latent  heat  fluxes  are  mostly  below  10 
Watts/m2.  Station  E  in  Region  II  is  a  notable  exception  where 
the  DMs  are  as  large  as  20  Watts/m2.  According  to  Table  X, 
Station  E  also  has  the  second  largest  mean  latent  3H  heat  flux 
at  124. 7  Watts/m2. 

The  difference  variances  appear  in  Figure  30.  Since  the 
DVs  at  all  averaging  periods  could  not  be  set  simultaneously  to 
0,  they  are  all  larger  than  those  that  appear  in  Appendix  D  for 
the  stresses  (the  heat  flux  DVs  are  indentically  0).  Up  to 
L=2.Q  days  the  DV  stress  error  for  all  components  is  bounded  by 
±  5. OX.  This  increases  to  ±  10. OX  for  the  majority  of  ships  by 
L=28.0  days.  An  exception  is  Station  D  L=28.0  days  which  will 
oe  investigated  later.  Op  to  L=7. 0  days  (with  the  exception  of 
Station  N) ,  the  DV  heat  flux  errors  are  all  bounded  between  ± 
10. OX.  From  L-14.0  to  1=28.0  days,  particularly  in  the  latent 
heat  flux,  the  CV  errors  become  guite  large  --  up  to  30. 7%  at 
Station  M.  These  are,  however,  all  substantially  lower  than  the 
60-80%  error  in  the  uncorrected  tests  outlined  in  Figures  7  and 
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sensible  heat  flux  latent  heat  flux 


Fig  ure  29  (a) 


Figure  29(b) 
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Fig-ure  29(c)  Linear  D.C. 


Figure  29. 
The  heat 
dPa.  D.C. 


The  difference  means  using  the  empirical  formula, 
flux  DMs  are  in  Watts/ra2  while  the  stress  D Ms  are  in 
denotes  drag  coefficient. 


21.  lhe  DV  error  in  the  stresses 
appears  to  oscillate.  Generally 
each  region  (ie.  L=. 25,  2.0,  4.0, 
lower  than  the  middle  points. 


and  sensible  heat  fluxes 
the  first  and  last  points  in 
and  28. C  days)  are  biased 
This  indicates  that  the 


adjustment  to  the  constant  (  )  terra  alone  of  Equation  5.  12  may 
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LATENT  HEAT  FLUX 


Figure  30  (a) 


Figure  30  (b) 
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Figure  30(c)  Linear  D.C. 

Figure  20.  The  difference  variances  using  the  empirical 
formula.  All  values  are  in  percent.  D.C.  denotes  drag 
coefficent. 


not  be  appropriate  and  that  an 
coefficients  may  be  necessary  to 
This  was  not  attempted  in  light  of 
uncorrected  DVs. 


adjustment  of  all  three 
force  all  the  DVs  to  zero, 
their  reductions  over  the 


The  residual  variances  appear  in  Figure  31 


For  nearly 
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every  ship,  averaging  period,  and  flux,  the  directly  applied 
empirical  formula  registers  an  improvement  over  the  RVs.  In 

the  stress  case,  the  improvements  are  marginal  --  being  limited 
to  about  1.05!.  Both  of  the  heat  fluxes,  however,  show  marked 
improvements.  At  L=28.0  days,  there  are  many  instances  where 
the  empirical  formula  residuals  are  about  one-half  cf  the 
values.  Ihis  is  quite  remarkable  when  one  considers  that  the 
formula  predicts  both  the  temporal  and  velocity  dependent 
corrections. 

Generally,  the  linear  drag  coefficients  RVs  increase  with 
averaging  period  in  Region  I  while  the  scatter  of  values  is  a 
fairly  constant  in  Region  II.  The  x  component  RVs  are  between 
6.0  -  12.05!  in  Region  II  while  -he  y  component  values  lie 
cetween  6.0  -  17.0?.  The  corresponding  constant  drag 

coefficient  values  are  about  one-half  of  the  linear  drag 
coefficient  values.  The  sensible  heat  fluxes  generally  follow 
the  pattern  of  the  stress  and,  with  the  exception  of  Station  N 
are  less  than  8.0?.  The  latent  heat  fluxes  are  somewhat  less 
accurate  with  all  stations  except  N,  E,  and  M  being  less  than 
18.0?.  Again  a  marked  reduction  from  the  uncorrected  values  has 
occurred. 

The  correlation  coefficients  appear  in  Figure  32.  Those 
points  having  low  RVs  have  high  correlation  coefficients  and 
vice-versa.  For  the  stresses  the  correlations  are 

marginally  larger  (ie.  better)  than  the  correlations,  having 

improvements  cf  several  parts  per  thousand.  More  substantial 
improvements  are  shown  in  rhe  heat  fluxes  where  improvements  of 
0.05  to  0.1  are  common. 
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Fig  ure  32  (a) 


Figure  32(b) 
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Figure  32(c)  Linear  D.C. 

Figure  32.  The  correlation  coefficients  with  the  ships' 
empirical  formula  applied.  D.C.  denotes  drag  coefficient. 


5.4  Geographical  Averaged  Results 

The  suns,  sums  of  squares,  and  the  sums  of  the  cross 
products  of  the  VA  and  3H  variates  were  saved  from  the  initial 
£  w<L  calculations  and  were  then  added  (excluding  Station  N)  to 
create  a  set  cf  tctal  statistics  for  all  the  ships.  The  non¬ 


linear  regression  (Equation  5.9)  was  then  applied  to  achieve  an 
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estimate  of  geographically  averaged  regression  coefficients 
(  ^  )  and  hence  a  single  set  of  formulae  applicable  to  most 
ship  locations.  The  results  appear  on  the  AVG  row  of  Appendix  H 
and  in  Table  XII. 

The  geographically  averaged  empirical  formula  (  ^  )  was 
then  used  to  determine  the  four  test  functions  at  each  location. 
The  complete  results  for  all  ships  including  Station  N  and  K, 
appear  in  Appendix  J.  Again  only  the  x  component  linear  drag 
coefficient  and  the  heat  fluxes  will  be  presented  . 

The  difference  means  are  shown  in  Figure  33.  In  general 
the  range  of  values  is  similar  to  those  of  the  results 
with  most  values  being  less  than  0.  1  dPa  for  the  stress,  and  5.0 
and  12.0  Watts/m2  for  the  sensible  and  latent  heat  fluxes.  The 
inapplicability  cf  at  Station  N  is  readily  evident  in  the 
abnormally  large  DM  values  for  all  three  fluxes  (see  Appendix 
J)  • 

The  difference  variances  appear  in  Figure  34  excluding 
Station  N  for  all  fluxes  and  Station  K  for  the  sensible  heat 
flux.  For  all  fluxes,  the  DV  values  have  a  range  approximately 
twice  that  of  the  £  e«  DVs.  The  oscillating  nehaviour  is  no 
longer  evident.  In  the  stress  case,  the  DVs  are  dispersed 
fairly  regularly  around  0.0;?.  The  heat  fluxes  ,  however,  have  a 
positive  bias.  Station  K  was  not  removed  from  averaging  in  the 
sensible  heat  fluxes.  The  effect  of  the  geographical  formula  on 
Station  K  is  reflected  in  its  abnormally  low  DV  values  in  the 
sensicle  heat  flux  (see  Appendix  J)  --  indicating  that  the 
geographical  averaged  formula  overestimated  the  3H  flux.  There 
are  two  possible  causes  for  the  positive  DV  bias.  The 
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TABLE  XII 

The  geographically  averaged  ships'  regression  coefficients  for 
the  form  =1*  <*(u2«- vZ)4/l  I  *  .  The  C  denotes  the 
constant  drag  coefficient,  the  L  denotes  the  linear  drag 
coefficient,  the  X  denotes  the  x  component,  the  Y  denotes  the  y 
component,  the  H  denotes  the  heat  fluxes,  the  S  denotes  the 
sensible  heat  flux  and  the  I  denotes  the  latent  heat  flux. 
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geographical  average,  ,  may  generally  reduce  the  residuals 
sufficiently  so  the  DV  values  begin  to  approach  the  RV  values; 
or  the  effect  of  including  Station  K  is  to  weight  the  statistics 
so  the  calculated  regression  coefficients  slightly  under¬ 
estimate  the  required  corrections.  Since  the  positive  bias  is 
also  felt  in  the  latent  heat  flux  where  Station  K  did  not  differ 
radically  from  the  other  locations,  the  former  is  the  probable 
cause  of  the  bias. 

Ihe  residual  variances  appear  in  Figure  35.  The 
geographically  averaged  values  are,  with  the  exception  of 
Station  N ,  remarkably  similar  to  the  individual  ship  empirical 
formula  7^  BV  values.  Fo-  the  stresses  and  the  sensible  heat 
fluxes,  the  RVs  at  all  averaging  periods  are  within  1.QX  of 
those  xisted  in  Appendices  D  and  G.  Station  K  sensible  heat 
flux  values  show  a  5.05?  increase  in  the  residuals  at  L=28  days 
over  the  individually  calculated  values.  As  noted  in  Chapter  IV 
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Figure  33.  The  difference  means  using  the  geographically 
averaged  empirical  formula.  The  heat  flux  EMs  are  in  Watts/m* 
while  the  stress  DMs  are  in  dPa.  DC  denotes  drag  coefficient. 
The  Station  N  values  are  omitted. 


,  the  Station  K  values  differ  quite  markedly  from  the  other 
snips'  values.  The  8. OS  residual  listed  in  Appendix  J  is  not, 
however,  markedly  larger  than  those  calculated  for  the  other 


ships. 


In  several  instances,  the  geographically  averaged  RVs  show 
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Figure  34.  The  difference  variances  using  the  geographically 
averaged  empirical  formula.  All  values  are  in  percent.  DC 
denotes  drag  coefficient.  Station  N  is  excluded  in  all  fluxes 
and  Station  K  is  excluded  from  the  sensible  heat  flux. 


an  improvement  over  the  individual  ship  RVs.  For  example,  at 
Station  D  L=  4. 0  days,  the  geographically  averaged  latent  heat 
flux  RV  is  12.2*  while  the  individual  RV  is  13.3*.  k  comparison 
between  the  respective  difference  variances  indicates  a  rise  to 
13.8*  from  -7.0*.  Thus  the  geographic  empirical  formula  has 
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Figure  35.  The  residual  variances  using  the  geographically 
averaged  empirical  formula.  All  values  are  in  percent.  Station 
N  is  excluded.  DC  denotes  drag  coefficient.. 


more  closely  approached  the  least  sguares  estimate  of  the 
Station  D  parameters  and  negated  the  DV  bias  correction. 

The  latent  heat  flux  geographic  >(  residuals  are 
sunstantially  larger  in  several  instances  than  the  individually 


calculated  values.  For  example,  at  Station  I,  the  geographic 
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residuals  at  L=23.0  days  are  27.6%  where  the  individual  formula, 
7^  ,  gives  18.0%.  In  all  other  cases,  the  geographic  residuals 
are  within  5.0%  of  the  individually  calculated  values. 

Ihe  correlation  coefficients  appear  in  Figure  36.  In  all 
cases,  the  geographically  averaged  correlations  are  about  0.02 
less  than  the  values. 

Ihe  geographically  averaged  empirical  formula  approach# 
with  two  exceptions,  appears  to  yield  nearly  identical  results 
to  tne  individual  ships  empirical  formula.  This  approach  is  not 
viable  for  any  of  the  fluxes  at  Station  N  and  may  not  be  the 
i/est  estimate  for  the  Station  K  sensible  heat  flux  correction. 
For  Station  D,  particularly  at  longer  averaging  periods,  the 
geographical  empirical  formula  shows  RV  improvement  over  the 
individual  ships  empirical  formula  because  no  attempt  is  made  to 
remove  the  DV  bias.  Ihis  is  also  reflected  in  the  general 
positive  DV  bias  indicating  that  the  geographical  average  may 
not  predict  the  long-term  3H  variance  quite  as  accurately  as  the 
individual  ships'  empirical  formula  (  ). 

5.5  Intrinsic  Temporal  Variations 

Ihe  latent  heat  flux  corrections  at  Station  N  are 
characterized  by  markedly  larger  residuals  and  lower  correlation 
coefficients  for  averaging  periods  in  Region  II.  When  the 
empirical  formula,  ,  is  applied,  at  L=28.0  days,  the 

residuals  are  as  high  as  58.9%  of  the  3H  variance.  Although  the 
sensible  heat  flux  residuals  are  larger  at  Station  N  they  are 
less  pronounced  than  fcr  the  latent  heat  flux. 

Ihe  required  transformations  may  be  inherently  non- 
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Figure  36.  The  correlation  coefficients  using  the 
geographically  averaged  empirical  formula.  DC  denotes  drag 
coefficient. 


stationary.  It  has  been  assumed  that  the  winter  corrections, 
for  example,  are  identical  to  the  summer  corrections.  Dorman 
(1574)  has  demonstrated  that  at  Station  N  there  are  large 
amounts  of  energy  at  the  annual  cycles  for  the  winds,  wind 


components 


sea 


surface  temperatures,  air  temperatures  and  dew 


Foint  tenperat ures.  Except  for  sea  surface  temperatures,  the 
character  of  all  component  variates  changes  aarkedly  from  winter 
to  suaaer.  Dcraan  also  showed  that  the  winter  clockwise  wind 
co»fcnent  power  is  lore  than  four  tines  that  ot  the  sunaer 
power.  Soae  of  the  scatter  that  is  evident  between  the  JH  and 
V A  variates  nay  be  due  to  seasonal  changes  in  the  regoired 
cor  rect  ions. 

One  aethod  to  determine  whether  the  annual  cycle  Flays  a 
dominant  rcle  in  the  required  corrections  is  to  examine  the  3H 
and  V A  and  residual  (3H  -  VA)  spectra.  If  a  significant  peak  in 
the  residual  spectra  occurs,  then  a  seasonally  dependent 
correction  aay  be  necessary.  To  this  end  Station  C  was  selected 
as  a  ccnttol  because  it  exhibits  particularly  stable  heat  flux 
residuals  at  L=28.0  days  --  7.6*  for  the  latent  heat  flax  and 
3.8*  for  the  sensible  heat  flux  as  coapared  to  58. 91  and  12.6X 
respectively  at  Station  N  after  the  ships  eapirical  focnula, 
V  ,  was  applied.  The  data  at  both  ships  were  averaged,  and  the 
lunar  aonthly  3H,  VA  and  residual  stresses  were  calculated.  The 
three  tiae  series  were  then  broken  into  9  two-year  blocks  at 
Station  C  and  10  two-year  blocks  at  Station  M. 

Each  block  of  data  was  detrended  and  a  cosine  taper  was 
applied  (following  Bendat  and  Piersol  (1971))  to  10.011  of  the 
data  at  each  end  of  the  blocks  to  reduce  ringing  fcon  the 
convolution  in  frequency  space  of  the  box  window  inherent  in  the 
data.  The  spectral  characteristics  of  the  box  window  and  the 
cosine  taper  can  be  seen  in  Bendat  and  P iersol  (197 1)  pp.  324- 
325.  The  cosine  taper  also  tends  to  reduce  the  power  of  the 
input  signal.  Consequently  the  variance  was  calculated  before 
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and  after  the  taper  and  all  the  output  spectra  multiplied  by  the 
appropriate  correction  at  each  frequency. 

The  twenty-six  points  within  each  block  were  then  fast 
fourier  transformed  by  the  (JBC  FOURT  programme  and  the  square  of 
the  magnitude  of  each  fourier  sine/cosine  pair  was  calculated. 
Ihe  spectral  density  estimates  for  all  blocks  were  added  and  the 
mean  power  density  estimate  was  found. 

The  error  in  the  mean  spectral  density  estimates  follows  a 

x. 

X-  distribution  with  2N  degrees  of  freedom  (Jenkins  and  Matts, 
1972)  where  N  is  the  number  of  blocks.  The  95%  confidence  zone 
is  determined  frcm: 

-  - 4  ^ 

%v*>,  o  <??>"  5,12 


where  V  is  the  number  of  degrees  of  freedom  (2N)  and 

XVe  lif  are  the  "X.  statistics  at  0.025  and  0.  975  probabilities 
and  9  degrees  or  freedom,  s*  is  the  spectral  density  estimate 
and  is  the  actual  unknown  spectral  density.  Equation  5.12 

is  correct  assuming  tnat  the  input  series  is  Uaussian  and  that 
tne  spectral  estimates  uetween  blocks  are  independent. 
Deterministic  signals,  ncn - s ta t iona r i t y ,  and  non-normality  of 
the  input  series  reduce  the  actual  degrees  of  freedom  causing  an 
increase  in  the  contidence  zones.  Sophisticated  techniques  nave 
neeii  developed  to  calculate  tne  effective  degrees  of  freedom. 
Inese  will  not  be  implemented  because  the  aim  is  to  determine 
oruy  if  outstanding  peak;;  t-xist  in  the  data.  Consequently  it 
wiri  re  assumed  that  is  ilso  the  <=ifectiv  -  degrees  of  freedom 
and  tne  factor  (v»  ,,  /*  1  will  r  -  iuoted.  For  *vo  legrees 

t.  0  s  ip  X-  ^  c  it>'  ^ 
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cf  freedoa  {ie.  one  independent  spectral  estimate) ,  this  factor 
is  approxiaately  (.25,4.0). 

The  spectra  of  the  uncorrected  latent  and  sensible  heat 
fluxes  at  Stations  C  and  N  appear  in  Figure  37.  In  all  cases, 
the  annual  cycle  in  the  3R  flux  is  evident.  At  Station  C,  the 
VA  pcwer  is  always  less  than  the  3H  power  and  if  one  considers 
the  3H  peaks  tc  a  he  a  factor  of  4.0  above  the  3H  noise  level, 
they  are  significant  to  less  than  2.0  degrees  of  freedoa. 

At  Station  N ,  however,  the  VA  power  densities  often  exceed 
the  3H  values  in  +he  latent  heat  flux,  consistent  with  the 
neqative  uncorrected  DV  values  shown  in  Figure  21.  In  the 
sensible  heat  flux,  the  VA  annual  cycle  is  well  within  the  VA 
noise  level  while  the  3H  annual  cycle  is  well  above  the  3H  noise 
level.  Furtheracre,  the  residual  series  power  density  for  both 
fluxes  at  the  annual  cycle  is  such  greater  than  the  residual 
noise  level.  In  fact,  the  annual  cycle  frequency  band  for  the 
latent  heat  flux  accounts  for  42%  of  the  total  power  of  the 
residuals.  The  processes  behind  this  power  nay  also  influence 
the  sewi-annual  cycle  causing  slightly  larger  power  densities  at 
this  frequency. 

The  individual  empirical  forwula  ,  >7  ,  was  applied  and  the 

( 

spectra  for  the  hear  fluxes  recalculated.  The  results  appear  is 
Figure  38. 

At  Station  c,  sensible  heat  flux,  the  VA  series  closely 
Batches  the  JH  series  d*  all  frequencies.  At  0.5  cycles/fear 
the  VA  series  sliqh*-ly  undei  {relicts  the  iH  series  while  at  3.5 
cycles/yeat,  it  slightly  over  predicts  the  1H  aeries.  The  over 
and  under  predictions  co u n t e i - ba 1  a  nee  accounting  for  the  lov 
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hi-jure  17.  ih-  ij«  ct  id  cf  the  uncor:ect?d  VA,  residual,  and  3H 
n<  <it  tlux  *  i  He  ■  enes  a*  Stations  c  and  N.  The  VA  series  is 
rriioatej  by  ♦  ,  t  he  ■  H  t  tr ies  by  x  and  the  residual  series  by 
^  •  I  he  t '  S  <of;  t  ;  ■!•=  i.ee  lieits  art  (.57,2.19)  for  Station  C  and 
til  i  or  s*  a*  i  ci  N. 
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Figure  38.  Ihe  spectra  of  the  corrected  heat  flux  for  the  3H , 
V A  and  residual  series  at  stations  c  and  N.  The  VA  series  is 
indicated  by  ♦  ,  the  3H  series  by  x,  and  the  residual  series  by 
O  •  ihe  confidence  limits  are  (.57,2.19)  at  Station  C  and 
(.  59,2.  08)  at  Station  N. 
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0.51  DV  level  shewn  in  Appendix  I.  The  residual  series  power  is 
quite  negligible  at  all  frequencies  compared  to  the  3h  series. 
The  latent  heat  flux  VA  series  also  closely  approximates  the  3U 
series  for  all  but  the  annual  cycle.  The  350  (Watt/a2) 2  years 
difference  at  this  frequency  accounts  for  a  large  portion  of  the 
3.01  difference  variance.  The  residual  series  shows  a  slight 
increase  at  the  annual  cycle  but  it  cannot  be  classified  as 
significant.  Thus,  even  though  there  is  substantial  power  in 
the  3H  and  VA  variates  at  the  annual  cycle  there  is  no 
indication  of  a  substantial  peak  in  the  residual  series. 

At  Station  N,  the  sensible  heat  flux  results  approximate 
the  3H  results  reasonably  well  at  all  freguencies  except  the 
annual  cycle  where  the  VA  results  underestimate  the  3H  results. 
At  8  other  harmonics,  the  VA  series  slightly  overestimates  the 
3h  series.  In  attempting  to  reduce  the  DV  error  to  0.0,  the 
variance  difference  at  the  annual  cycle  may  have  been  increased 
to  the  detriment  of  the  other  frequencies.  The  residual  series 
shows  a  slight  rise  at  the  annual  cycle  but  is  probably  not 
significant  in  comparison  to  the  error  at  the  other  freguencies. 
In  the  latent  heat  flux  at  Station  N ,  the  VA  series 
underestimates  the  3H  spectral  densities  at  most  frequencies. 
At  the  annual  cycle,  the  residual  power  density  is  more  than  4 
times  that  at  any  other  frequency.  Thus  seasonal  corrections 
■ay  be  required  at  Station  N  in  the  latent  heat  flux  and  this 
»ay  in  part  account  for  its  larger  residual  variances.  The 
residual  power  density  is  approximately  equal  to  the  sum  of  the 
VA  and  3H  powers.  It  is  demonstrated  in  Appendix  A. 4  that  this 
can  occur  when  a  there  is  a  90°  phase  shift  between  the  JH  and 
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VA  variates. 

Ti,e  time  series  of  the  residual  latent  heat  flux  at  Station 
N  with  yj  applied  is  shown  in  Figure  3S(a).  The  annual  cycle 
is  evident.  In  contrast,  the  Station  C  residuals  are  in  Figure 
39(b).  As  suggested  in  the  spectra,  there  are  no  dominant 
signals. 

An  inspection  of  the  data  at  Station  N ,  L=23.0  days,  latent 
heat  flux  revealed  that  grouping  months  1-5  and  11-13  inclusive 
reflected  a  positive  bias  in  the  residuals  while  months  6-10 
inclusive  reflected  a  negative  bias  in  the  residuals.  The  data 
were  then  arranged  in  these  groups  and  a  winter  and  summer 
empirical  formula  calculated  for  L=28.0  days.  The  results 
appear  in  Table  XIII.  The  summer  values  of  about  -1.8  are 
substantially  lower  (more  negative)  than  any  other  values.  The 
latent  heat  residuals  were  calculated  using  these  values  and 
were  reduced  to  31.9?  and  the  sensible  heat  residuals  were 
reduced  to  6.0?. 

The  spec*ta  for  the  seasonally  corrected  fluxes  at  Station 
N  appear  it  r : a  it-  «o.  3H  an)  VA  series  VA  series  now  match 
more  closely  at  t  t.  : :  r  i  :  1  y  le.  A  comparison  with  Figure  38 
reveals  ‘ha*  i  i. .  .  a  . a  the  power  of  the  residuals  at 
2.0  cyci-v./.  !  .  ;at  reordering  of  the  relative 
positrons  :  f.  •  :--''Tal  densities  has  occurred  at 
most  ttvjuc'ii  :  .  -  .  •  <t  i  a  i  density  of  the  residual  series 
at  the  annual  cvri-  :  s  well  within  the  scatter  of  the  other 
fregaencies.  The  j  r  of  the  residuals  at  all  freguencies 
other  than  tne  annual  cycle  lias  remained,  however,  virtually 
unchany  ed. 
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TABLE  XIII 

The  summer  and  winter  formula  coefficients  required  at  Station  N 
for  the  heat  fluxes  at  L=2S.O  days. 

i - t - t - - - — i 

J  TYPE  |  WINTER  i  SUMMER  I 

V - + - + - + - + - ^ 

i  i  oi  i  /S  i  i  rs  i 

j. - - + - .( - (. - j 

1  Sensible |  5.337  |  -1.368  |  4.399  |  -1.730  | 

|  Latent  |  4.997  |  -  1.  373  1  3.  676  |  -  1.824  | 

! _ X _ J _ X - L - 1 

Station  D  has  the  largest  residuals  for  the  y  component, 
linear  drag  coefficient  and  the  second  largest  residuals  for  the 
x  component  linear  drag  coefficient.  Consequently,  the  spectra 
at  this  station  were  examined  at  L=28.Q  days.  Figures  41(a) 
and  (n)  are  the  x  and  y  ccmfcnent  spectra  calculated  from  the 
uncorrected  time  series.  The  x  component  3h  and  residual 
spectra  both  show  distinct  increases  in  spectral  density  at  the 
annual  cycle  while  the  y  component  is  guite  white.  In  both 
components  the  residual  spectra  are  larger  at  all  frequencies 
than  the  VA  spectra  indicating  that  the  error  is  greater  than 
the  actual  VA  estimate  of  the  power.  This  indicates  that  severe 
errors  may  arise  if  geostrophic  winds  from  barometric  pressure 
maps  averaged  over  periods  cf  greater  than  perhaps  four  days  are 
used  to  calculate  stress  spectra.  One  expects  much  smaller 
errors  in  the  spectra  calculated  by  Willebrand  (  1 978)  since  the 
errors  for  the  28  day  averaging  period  shown  here  is  much  larger 
than  those  for  the  12  hourly  averaging  period  upon  which  he 
based  his  study. 

Figures  41(c)  and  (d)  show  the  VA  and  3H  spectra  with  ^ 
applied  to  the  VA  series.  In  the  x  component  there  is  a 
definite  mis-match  between  the  3H  and  VA  spectra  at  the  atnual 
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cycle  while  the  residual  spectrum  shows  a  slight  peak;  At  all 
other  frequencies,  hcwever,  the  3H  and  VA  spectra  match  quite 
well.  After  correction,  the  y  component  3H  and  VA  match  fairly 
well  at  all  frequencies  while  the  residual  series  is  much  lower 
in  power.  Since  the  y  component  power  scale  is  1/2  of  the  x 
component  power  scale,,  the  variance  of  the  residuals  of  the  y 
scale  are  somewhat  less  than  1/2  of  those  of  the  x  component. 
Taole  V  indicates,  however,  that  the  actual  3H  variance  in  the  y 
component  at  Station  D  is  only  1/10  of  that  in  the  x  component. 
Thus  the  generally  larger  relative  residual  variances  shown  in 
the  y  component  for  the  linear  drag  coefficient  may  be  partially 
attributed  to  dividing  by  the  smaller  3H  variance. 

The  geographical  empirical  formula,  7[  ,  was  applied  to 
the  linear  drag  coefficient  stress  at  station  D  and  the 
resulting  spectra  appear  in  Figure  42.  The  residual  spectra  are 
virtually  identical  to  those  obtained  when  the  individual 
correction,  ,  was  applied.  The  VA  spectra,  however, 
systematically  underestimates  the  3H  series  spectral  densities 
at  all  frequencies.  This  is  consistent  with  the  increased 
difference  variances  induced  by  application  of  the  geographical 
empirical  formula. 

The  Station  N  residual  series,  latent  flux  shows  a  large 
amount  of  power  at  the  annual  cycle  which  can  be  reduced  by 
applying  a  semi-annual  correction.  This  suggests  that  a  time 
dependent  correction  may  be  required  for  the  stresses  at  Station 


N.  The  individual  Station  N  stress  corrections  for  the  linear 
drag  coefficient  were  applied  and  the  results  are  shown  in 


1 
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density]  has  been  increased  by  an  order  of  magnitude  over  that 
for  Station  E  consistent  with  the  lower  variances  reported  in 
ladle  V.  There  is  a  slight  mis-match  of  power  between  the  3H 
and  VA  series  as  well  as  a  slight  rise  in  the  y  component  stress 
power  densities  at  the  annual  cycle.  The  increase  in  the 
residual  power  is  not  significant.  Consequently  a  monthly  or 
seasonal  dependence  in  the  corrections  may  not  be  required  for 
the  stresses. 
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CHAPTER  VI  MIND  ROSE  HEASUREHENTS 

6# 1  Introduction 

Over  much  of  the  Worlds's  oceans,  the  only  source  of 
climatic  intonation  is  from  weather  logs  of  ships  of 
opportunity.  The  Office  of  the  Chief  of  Naval  Operations,  US 
Navy  has  compiled  and  smoothed,  in  the  foraat  of  vind  roses,  the 
aetecroloqical  observations  of  weatherstations  and  ships-of- 
cppcrtunity  for  each  of  the  twelve  months  of  the  year  as  well  as 
a  ccnglc«erate  wind  rose  tased  on  all  the  available  data. 

The  wind  rose  is  constructed  by  sorting  the  data  into  eight 
directions  and  four  Beaufort  categories.  They  are  then 
presented  as  the  relative  frequency  of  wind  events  within  the 
four  direction/Beaufort  groups  and  as  the  relative  frequency 
within  nine  Beaufort  classes  irrespective  of  direction.  Further 
details  can  be  fcund  fro«  the  Marine  Climatic  Atlas(HCA). 

He  1  lerman ( 1965  and  1967)  has  used  wind  roses  to  estiaate 
world  cliBatic  wind  stress  patterns.  Using  a  discrete  and  a 
continuous  stress  calculation  technique  and  three  different  drag 
coefficients  to  coipute  stress  values,  he  estimated  the  wind 
stress  curl  and  calculated  different  values  of  the  wind  driven 
western  boundary  current.  He  concluded  that  the  selection  cf 
drag  coefficient  had  a  such  larger  influence  on  the  magnitude  of 
the  western  boundary  current  than  errors  arising  from  the  type 
of  calculation  of  stress  from  wind  roses. 

The  weatherstation  data  offered  an  opportunity  to 
investigate  in  detail  the  inaccuracies  arising  in  calculation  of 
wind  stress  frc«  wind  roses. 


To  facilitate  discussion,  the  nomenclature  of 
He 1 1€ nan  ( 196 5)  will  be  used.  The  variable  t  .  mill  refer  to 

O 

the  leiative  frequency  of  wind  events  in  Beaufort  category  i  and 
direction  j  while  Pi  will  refer  to  the  relative  frequency  of 
wind  events  in  Beaufort  cateqory  i  irrespective  of  direction. 
Note  that  Helleman  labelled  his  Beaufort  categories  fron  i  =  0 
-  11  where  the  present  study  will  use  intervals  i  =  1  12. 
Beaufort  category  13  (>  33.5  i/sec)  will  be  arbitrarily 
incorporated  into  class  12. 

The  original  wind  direction  was  digitized  in  10°  increments 
giving  36  directions.  The  Marine  Climatic  Atlas  was  formatted 
into  8  directions.  Because  eight  is  not  a  factor  of  36,  the 
intervals  between  the  two  wind  roses  do  not  coincide  and  a  bias 


occurs  in  the  constructions  of  the  8  direction  wind  rose  as 


outlined  by  Lea  and  Hel vey  ( 1 97 1)  . 

To  compensate  for  this  bias,  a  linearly  piecewise 
continuous  36  direction  distribution  for  each  Beaufort  category 
was  constructed.  Hellerman  used  a  sinilar  technique  to 
interpolate  wind  speeds;  we  adapted  the  technique  to  interpolate 
directional  distributions.  For  each  of  the  36  direction 


intervals,  f  was  assumed  to  be  centred  on  the  j  direction  (eg. 
5.0«).  The  relative  frequency  was  divided  by  the  width  of  the 
direction  category  (in  this  case  10®)  to  arrive  at  a  total 
frequency  density  for  each  direction/Beaufort  category.  The 


frequency  densities  at  the  ends  of  each  category 
determined  by  interpolation  between  adjacent  centre 


«*re 
points  as 


shown  in  Figure  49.  The  centre  points  were  then  readjusted  to 
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total  JfcO0  was  then 

divide  d 

in’-o  mocks 

of 

4  S 

0 

w  i  *  h 

*  he 

interval  centred  on  0.0°  and  the  area  hit  Inn  each  bloc 
calculated  to  determine  the  eight  direction  frequency  t  i:.  t 
The  MCA  format  wind  rose  groups  ’he  Beaufort  interval 
follows  : 

•  3  J  *  "i  •  *  7  v  *  j  I  f  c->  * 

^1,1.  '^3  •  ^  f  •  f  t  •  ^1  •  5  F|#  * 


1:  i=1,e  and: 


6.4 


Ignoring  the  intrinsic  error  of  the  transfer  coefficients 
and  accepting  those  used  in  the  previous  chapters,  three 
immediate  r  trees  of  errcr  can  be  identified  as  follows: 

a.  err -rs  in  the  directional  construction  of  the  wind  rose; 

b.  errors  in  wind  magnitude  grouping;  and 

c.  differences  in  calcuiation  techniques. 

Assuming  that  the  36  direction  wind  comprises  the  initial 
data  set,  the  reduction  to  8  directions  induces  errors  to  tne 
stress  calculations  in  twc  ways.  First,  since  the  linear 
piecewise  continuous  d ist ri out icn  is  an  approximation  of  the 
actual  continuous  distribution,  the  assumed  overlap  calculations 
may  not  be  accurate.  Second,  the  reduction  from  36  to  8 
directions  reflects  a  reduction  in  the  numoer  of  degrees  of 
freedom  of  the  distribution  (and  nence  th<=  total  information 
available  about  the  statistical  structure  of  the  system)  oy  a 
factor  greater  than  four.  A  wind  spee 5  averaging  is  inherent  in 
the  reduction  cf  number  of  directions  and  should  result  in  a 
decrease  in  accuracy  cf  the  stress  c a lc u la* ions . 

Errors  due  to  wind  speed  smoothing  arise  in  a  similar 


manner  to 

the  errors 

due  tc  di 

rectron  smoothing. 

ihe  original 

data  tor 

Stations  A 

through 

N 

are 

digitized  to 

1.0  knot  (0.513 

a/sec)  . 

Ihe  initial 

grouping 

to 

1 1 

Beaufort  intervals  markedly 

reduces  t 

he  number  of 

deg  rtes 

of 

f  re 

eaom  of  the 

distrioution  and 

ltd uces 

the  vector 

averagin 

of 

t  h  e  wind. 

The  MCA  format 

increases  the  wind  smoothing  oy  a  t irther  reduction  in  grouping 
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from  12  Beaufort  categories  per  direction  tc  four.  The  Modified 
Scripps  technique  attempts  to  overcome  this  reduction  by 
expansion  cf  the  MCA  format  to  12  Beaufort  categories  whereas 
riellerman  expands  through  interpolation  of  the  given  frequency 
densities. 

The  Modified  Scripps  calculation  is  essentially  a  discrete 
calculation  in  which  it  is  assumed  that  the  frequency  densities 
witnin  each  Beaufort  cateqory  should  have  little  effect  on  the 
final  flux  results.  The  Hellerman  tecnnique  assumes  that  the 
frequency  densities  are  continuous  and  that  an  attempt  to 
interpolate  the  variations  within  each  Beaufort  category  may  be 
more  erfective.  A  discrete  and  continuous  flux  calculation 
technique  can  be  applied  with  facility  at  all  levels  of  wind 
rose  construction. 

The  total  data  sample  was  divided  into  actual  monthly 
groupings  and  the  stresses  and  neat  fluxes  were  directly 
calculated.  Tc  simulate  irregular  samplings  by  ships-of- 
opportunity,  any  3H  reading  which  contained  erroneous  data  for 
any  variacle  was  rejected  from  calculations.  Wind  roses  were 
then  constructed  in  three  stages  including  monthly  and  total 
wind  groupings.  The  36  direction,  12  Beaufort  interval  was 
initially  constructed.  The  stress  was  then  calculated  by  the 
discrete  and  ccntinucus  methods  to  determine  the  loss  in 
accuracy  from  wind  speed  smoothing  alone.  An  8  direction,  12 
Beaufort  interval  wind  rose  was  then  constructed  from  the 
previous  wind  rose  and  the  stress  calculated  d y  both  techniques 
to  indicate  the  loss  in  accuracy  due  to  direction  smoothing, 
iii. ally  a  MCA  type  wind  rose  of  8  directions  and  4  Beaufort 


135 


intervals  was  constructed  and  the  stress  calculated  by  the 
Hellerman  continuous  technique  and  the  Modified  Scripps 
procedure  to  determine  the  loss  of  accuracy  due  to  further  wind 
speed  smoothing  and  the  accuracy  one  can  expect  from  Marine 
Climatic  Atlas  wind  roses. 


In  order  to 

facilitate 

identification. 

the 

type 

of 

calculation  will 

be  identified 

by 

the  number  of 

direc 

tions 

and 

number  of  Beaufort 

intervals  of 

the 

wind  rose  and 

the 

type 

of 

flux  calculation. 

The  type 

of 

flux  calculation 

will 

be 

abbreviated  as  :  (1)  D  -  discrete,  (2)  C  -  continuous  (3)  MSc 
modified  Scripps.  Thus  the  36-12-D  stress  will  be  the  stress 
calculated  from  the  36  direction  12  Beaufort  interval  wind  rose 
by  the  discrete  technique.  The  discrete  technique  applies 
Equation  6.2  directly.  The  MSc  method  involves  the  expansion  of 
the  four  Beaufort  categories  to  12,  then  the  application  of 
Eguation  6.2.  The  abbreviations  and  symbols  to  be  used  in  the 
stages  of  calculation  of  the  wind  rcses  appear  in  Table  XIV. 

Tne  stress  differences  were  quantified  by  determining  the 
absolute  differences  in  stress  magnitude  and  direction  between 
the  wind  rose  stress  and  the  3H  stress.  The  heat  fluxes  were 
determined  by  calculating  the  monthly  average  air-sea 
temperature  and  humidity  differences  directly.  The  wind  roses 
were  used  to  calculate  the  monthly  averaged  wind  speeds.  The 
wind  rose  heat  flux  was  then  the  product  of  the  wind  rose  wind 
speed  and  the  monthly  average  temperature  or  humidity 
difference.  The  difference  between  the  directly  determined  and 
wind  Lose  heat  fluxes  then  quantified  the  inaccuracies. 

The  differences  were  determined  for  each  monthly  group  and 
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TABLE  XIV 

Coding  for  tne  wind  rose  levels  calculated  in  Chapter  VI. 


r  - - 

i  DIIiZlTIONS 

i 

i  _ _ 

T - - 

|  SPEEDS 

1 

-_L  _ 

-  r - 

1  CALCULATION 
ITECHNIUUE 

4  . 

t  — - - 

| ABBREVIATION 

1 

~T“ 

1 

1 

SYMBOL 

— i 

1 

1 

i 

36 

t 

1 

12 

(  Continuous 

| 36-  12-C 

1 

1 

i 

36 

1 

12 

| Discrete 

| 3  6-  12-D 

1 

© 

i 

i 

8 

1 

12 

I  Con t in  uou  s 

| 6-12-C 

1 

X 

i 

i 

8 

1 

12 

| Discrete 

| 8-  12-D 

1 

A 

i 

i 

8 

1 

MCA 

|Cont inuou  s 

I  MCA 

1 

i 

i 

i 

i 

8 

1 

1 

1 

MCA 

I  Icontinuous 

| M  Sc  |  MC  A  MSC 

1  i 

1 

1 

* 

i 

i 

1 

L - - 

— 

1 

1 

1 

for  the 

composite 

record 

wind  rose  for 

each  ship. 

The  monthly 

differences 

were 

then 

averaged  in 

two  manners 

• 

First , 

the 

annual 

cycle 

for 

each 

ship  was  averaged  and 

the 

standard 

deviations  of  the  monthly  discrepancies  determined.  For 
example,  the  discrepancies  from  January  through  December  for 
Station  A  were  averaged  tc  determine  a  typical  error  for  Station 
A.  Inis  averaging  will  be  referred  to  as  the  ship  average  error 
and  the  corresponding  standard  deviation,  the  ship  error 
standard  deviation.  The  second  averaging  was  performed  between 
ships  for  each  month.  For  example,  the  January  errors  only  for 
all  ships  were  summed,  averaged  and  the  standard  deviations 
determined.  This  technique  should  indicate  if  the  errors  are 
induced  by  the  annual  cycle  and  will  be  referred  to  as  the 
monthly  average  and  monthly  error  standard  deviation. 


6.3  ties  ult s 

The  stress  and  heat  flux  components  were  calculated  through 
the  six  stages  previously  outlined.  In  all  cases  the  errors 
involved  were  generally  much  smaller  than  the  difference  mean 
errors  discussed  in  Chapters  III,  IV,  and  V. 


In  most  cases  the 
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error  in  stress  magnitudes  were  less  than  ±  0.05  dPa;  the  errors 
in  direction  were  generally  less  than  5.0°  ;  and  the  heat  flux 
errors  were  generally  bounded  between  ±5.0  Watts/a*. 

Stress  Magnitude  Errors 

Figures  45  (a)  and  (b)  contain  the  ship  mean  and  error 
standard  deviations  for  the  constant  drag  coefficient 
respectively.  Fiaures  45  (c)  and  (d)  contain  the  monthly  aean 
and  error  standard  deviations  respectively  for  the  linear  drag 
coefficient. 

In  all  cases  ♦  he  36-12-d  calculation  is  biased  slightly 
neqatively  at  about  0.01  dPa  indicating  that  this  technique 
consistently  eve resti nates  the  actual  stress.  This  is 
accompanied  by  a  standard  deviation  of  siarlar  order  of 
■aqnitude  suggesting  that  the  bias  is  statistically  significant. 
It  is  doubtful,  however,  •■hat  errors  of  this  aagaitude 
contribute  greatly  to  induce  errors  in  oceanic  circulation 
calculations. 

In  both  the  ship  and  monthly  means,  the  36-12~C 
calculations  are  virtually  0.0  dPa.  The  error  standard 
deviations  are  roughly  eguivalent  to  the  36-12-D  technique. 

The  8-12-C  has  a  mean  error  of  about  0.02  dPa  larger  than 
the  8-12-D  calculations  with  nearly  identical  error  standard 
deviations.  Furthermore,  both  cases  show  a  consistent  reduction 
of  C.01  -  0.02  dPa  from  the  respective  36  direction  calculation. 
Thus  an  effect  c£  reduction  of  direction  degrees  of  freedom  is 
to  induce  a  slight  underestimation  of  the  stress.  It  is 
interesting  that  this  induced  underestimation  in  the  discrete 


IFFERENCE  IN  DPfl^  CONSTANT  MAG  DIFFERENCE  IN  DPR 
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calculation  re-adjusts  a  systematic  overestimation  for  36 
directions  to  an  over-all  mean  of  nearly  0.0  dPa  in  the  8 
direction  case.  In  the  continuous  calculation,  the  means  are 
consistently  underestimated  by  about  0.02  -  0.03  dPa. 

The  MCA  continuous  and  discrete  calculations  lead  to 
curious  results.  The  continuous  calculation  mean  errors  (with 
several  exceptions) ,  appear  to  underestimate  the  total  stress  to 
about  the  same  degree  as  the  8-12-C  technique.  The  HC. A  MSc  mean 
values  are,  however,  lower  than  the  8-12-D  technique.  Indeed, 
during  the  summer  months  the  MCA  MSc  technique  overestimates  the 
"actual"  stress  to  an  even  greater  degree  than  the  36-12-D 
technique.  The  standard  deviations  for  the  discrete  and 
continuous  calculations  are  nearly  identical  and  display 
distinct  increases  (up  to  0. 1  dPa)  in  the  winter. 

The  HCA  format  does  not  resolve  the  gale  force  winds  into 
direction  categories.  To  achieve  this,  Hellerman  divides  1/6  of 
the  total  gale  force  winds  into  Beaufort  categories  11  and  12 
irrespective  of  season  and  location.  Since  the  stress  is 
calculated  as  the  square  of  the  wind  speed  (a  cubic  io  the 
linear  case) ,  even  small  gale  force  wind  frequency  densities  can 
make  disproporticna tely  large  contributions  to  the  total  stress. 
Thus,  if  the  actual  ratio  of  Beaufort  categories  11  and  12  is 
substantially  less  than  1/6  then  a  fictitiously  high  stress  will 
result.  For  example,  at  Station  J  where  the  discrete 
calculation  yields  a  particularly  high  value  in  January,  the 


qale 

force  winds  account 

for 

18.0  % 

of  the  distrib 

ution. 

The 

ratio 

of  Beaufort  categories 

1 1 

and  12 

to  the  total 

gale 

force 

winds 

is  only  11.3%  and 

not 

16. 7%. 

Furthermore, 

in  June  and 
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July  there  are  no  recorded  values  in  Beaufort  categories  II  and 

12. 

All  errors  calculated  by  all  the  techniques  a r<  v  y  snail. 
The  reduction  xn  directional  degrees  of  freedoa  leads  to  a 
slight  increase  in  the  underestimation  of  the  "true"  stress. 
The  discrete  calculation  consistently  yields  larger  stress 
magnitude  estimates  than  the  continuous  calculation  for 
equivalent  degrees  of  freedom.  Inaccuracies  associated  with  the 
KC  A  techniques  nay  be  due  to  inadequate  interpolation  of  the 
gale  force  winds. 

Stress  Direction  Errors 

The  ship  linear  angle  difference  neans  and  standard 
deviations  appear  in  Figures  46  (a)  and  (b)  respectively  while 
the  monthly  constant  angle  difference  neans  appear  in  Figures  46 
(c)  and  (d)  . 

The  mean  angle  differences  are  all  less  than  5.0®.  The 
standard  deviations,  with  the  exception  of  Stations  D  and  K, 
linear  draq  coefficient,  are  also  less  than  5.0®.  The  neans  for 
the  constant  drag  coefficient  as  a  function  of  month  are  biased 
positively  while  the  lineal  drag  coefficient  biases  appear  less 
evident. 

The  mean  errors  appear  to  be  dependent  upon  degrees  of 
freedom  and  not  upon  discrete  or  continuous  calculation 
techniques.  The  36-  12-C,  36-12-D,  8-12-C,  and  8-12-D  techniques 
are  virtually  identical  in  both  neans  and  standard  deviations. 
The  two  MCA  calculations  are  also  very  close  in  their  nean 
values  although  the  parallel  is  such  less  evident  in  the 
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figure  46(c)  Constant  CC  Figure  46(d)  Constant  DC 

Figure  46.  The  direction  mean  errors  and  error  standard 
deviations  between  the  actual  stress  and  wind  rose  stress  as  a 
function  cf  ship  and  month.  The  T  is  the  composite  wind  rose 
based  on  all  the  data.  DC  denotes  drag  coefficient. 
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standard  deviations.  Where  the  two  MCA  standard  deviations 
diverge,  the  continuous  calculation  shows  lore  fluctuations. 

Heat  Fluxes 

The  errors  in  the  sensible  heat  fluxes  appear  in  Figure  47. 
Errors  are  very  snail  --  the  leans  are  generally  bounded  by  ± 
2.0  Watts/m*.  Within  each  ship  ct  month  category  the  standard 
deviations  for  all  calculations  techniques  are  virtually 
identical.  The  reduction  from  36  to  8  directions  appears  to 
have  a  negligible  effect  on  the  mean  error  for  either  the 
discrete  or  continuous  calculations.  There  is,  however,  a 
slight  positive  bias  (or  underestimation)  in  the  mean  error  for 
the  36-12-C,  36-12-L.,  8-12-C,  and  3-12-D  calculations. 

Wind  speed  grouping  according  to  the  MCA  foraat  induces  a 
small  mean  error.  The  continuous  MCA  calculation  consistently 
underestimates  he  sensible  heat  flux  in  comparison  to  the  other 
calculations  while  the  MCA  HSc  calculations  yield  a  consistent 
overest imation  in  relatLon  to  the  other  techniques.  The  effect 
of  this  slight  o verestimat ion  is  that  the  MCA  HSc  calculation  on 
the  average  has  zero  bias  in  the  means  for  both  the  ship  and 
monthly  average  calculations. 

The  latent  heat  flux  results  appear  in  Figure  48.  Bany  of 
the  observations  of  the  sensible  heat  flux  ace  also  applicable 
to  the  latent  heat  flux.  There  is  little  difference  in  the 
means  between  any  calculation  for  the  direction  reduction  and 
the  major  discrepancies  occur  for  the  MCA  foraat.  The 
continuous  MCA  calculation  again  underestimates  the  latent  heat 
flux  in  relation  to  the  other  techniques  while  the  BCA  BSc 
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Figure  48.  The  mean  errors  and  error  standard  deviations 
between  the  actual  latent  heat  fluxes  and  the  wind  rose  latent 
heat  fluxes  as  a  function  of  ship  and  month.  The  T  refers  to 
the  conglcaierate  wind  rose. 
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calculation  again  overestiaates  the  latent  heat  flux.  The 
standard  deviations  for  all  technigues  are  virtually  identical. 

The  nean  latent  heat  errors  are  largely  bounded  by  t  5 
Hatts/a2  with  a  general  negative  bias  for  for  all  technigues  in 
both  the  aonth  and  ship  lean  errors.  The  effect  of  the  negative 
bias  is  that  the  flC  A  continuous  calculation  now  provides  on 
average,  the  sost  accurate  estimates  of  the  latent  heat  flux. 
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CH4PTE8  VU  SUNHAh*  *12  CONCIJSUNS 

Time  secies  of  varying  lenqths  from  nine  weatherships  in 
the  North  Atlantic  0c»  an  and  fro*  two  weatherships  on  the  North 
Pacific  Ocean  were  examined.  The  data  consisted  of  three-hoiniy 
readings  of  wind  component  velocities,  barometric  pressure,  air 
temperature,  dew  point  temperature,  and  sea  surrace  temperature 
from  which  were  calculated  three-hourly  air  densities,  air-sea 
temperature  differences,  air-sea  humidity  differences,  x  ana  y 
component  wind  stresses  {using  both  constant  and  linear  drag 
coef ficients) ,  sensible  and  latent  neat  fluxes.  The  calculated 
quantities  and  the  wind  component  velocities  were  averaged  oter 
varying  periods.  The  fluxes  were  then  recalculated  using  the 
appropriate  averaged  constituent  variates  to  arrive  at  a  vector 
averaged  ( V  A )  estimate.  The  VA  flux  estimates  were  then 
compared  to  the  directly  averaged  (called  3H)  fluxes. 

A  systematic  underprodiction  of  the  3H  air-sea  fluxes 
occurs  when  they  are  estimated  through  VA  parameters.  The 
required  transformations  to  correct  for  these  losses  are  similar 
for  all  fluxes  and  are  dependent  upon  the  ? A  wind  speed  and  upon 
averaging  period.  (Chapters  III  and  IV) 

Four  test  functions:  the  difference  means,  difference 
variances,  residual  variances,  and  correlation  coefficients  were 
defined  in  Fquations  2.16  -  2.19  to  quantify  the  differences 
between  the  VA  and  3H  fluxes.  In  all  cases  the  difference 
means,  difference  variances,  and  residual  variances  increased 
roughly  exponentially  with  averaging  period  while  the 
correlation  coefficients  showed  a  similar  decrease.  The 


residual  variance  measured  the  point-by-point  differences  and 
had  maximum  unccrrected  values  at  L=28.0  days  of  30. 0-40. OX  for 
the  constant  drag  coefficient,  40.0-tjO.JX  for  the  linear  drag 
coefficient,  25. 0-45. OX  for  the  sensible  heat  flux  and  40. QX* 
tor  ♦he  lament  heat  flux.  (Sections  3.1  ana  4.2) 

A  simple  analysis  established  a  lower  bound  of  the  point- 
fcy-pcin*  magnitude  reductions  in  the  constant  drag  coefficient 
stress,  Sj(L).  They  were  proportional  to  the  variance  of  the 
wind  speed  lost  through  averaging,  and  inversely  proportional  to 
the  calculated  vector  averaged  wind  speed  squared.  Although  not 
exact,  *-he  simple  model  aided  in  explaining  many  distributional 
features  of  the  3H-VA  discrepancies.  The  fact  that  the  point- 
by-point  direction  differences,  8  j  ( L)  ,  averaged  to  0.0®  for 
moderate  wind  speeds  indicated  that  a  single  transf ornation 
would  suffice  for  both  the  x  and  y  stresses.  (Section  J .3) 

The  most  precise  method  of  correcting  for  the  3H-VA 
discrepancies  was  to  linearly  regress  the  two  variates  on  a 
Peaufort  cateqory/averaginq  period  basis  arriving  at  a 

y 

transformation  .  flaried  improvements  were  found  in  all  four 
test  quantities.  The  residuals  increased  exponentially  to  3.0 
days  and  then  remained  fairly  constant  from  3.0  to  28.0  days. 
At  L=28.0  days,  the  constant  drag  coefficient  residuals  were 
reduced  to  5. 0-8. OX,  the  linear  drag  coefficient  residuals  were 
reduced  to  10. 0-15. OX,  the  sensible  heat  flux  residuals  were 
reduced  to  5.0-10.0%  and  the  latent  heat  flux  residuals  were 
reduced  to  10. 0-40. 0*%.  Significant  (>10.0X)  t ransf or mations 
were  required  for  all  wind  speed  categories  for  averaging 
periods  greater  than  0.5  days.  (Sections  3.7  and  4.5) 
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Due  to  aathemat ica 1  aspects  of  the  regression,  it  was  found 
that  a  positive  bias  in  the  difference  variances,  exactly  equal 
to  the  residual  variances  was  induced.  Because  correlations 
were  generally  high  between  the  3H  and  VA  variates,  this  bias 
could  be  corrected  without  markedly  increasing  the  residuals. 
(Sections  3.5  and  4.3) 

Plots  of  indicated  that  the  transformations  may  be 

quite  reqular  ever  varying  averaging  periods  and  varying  wind 
speeds  for  all  fluxes.  A  non-linear  rejression  of  the  form 
X.  =1  ♦  ^[\i2*v2)^2L  was  performed  and  an  analysis  of  the 

test  quantities  revealed  that  application  of  this  formula  to  the 
VA  fluxes  produced  identical  results  to  the  £ kf  test  results. 
The  best  method  of  application  of  the  formula  was  to  allow  it  to 
correct  each  VA  estimate  on  a  point-b y- point  basis.  (Section 
5.2) 

The  similarity  of  the  results  for  all  fluxes  indicated  that 
the  vector  averaging  of  the  wind  in  general,  and  the  loss  of 
wind  variance  information  in  particular,  determined  the  extent 
of  the  reduction  for  all  fluxes.  This  was  confirmed  at  Station 
N  where  a  10  fold  reduction  in  wind  speed  variance  accounted  for 
markedly  reduce!  transformations  requited  tor  both  drag 
coefficient  ca leu  la t i ons ,  and  for  Doth  heat  fluxes.  The 
discrepancy  at  Station  K,  sensible  heat  flux,  remains  an  enigma. 
(Sections  3.1,  3.4,  and  4.4) 

It  was  shewn  by  direct  observation  that  all  ships  excluding 
Station  N  for  all  fluxts  and  Station  K  for  the  sensible  heat 
flux  required  similar  transfor ma t ions.  The  directly  calculated 
difference  means,  residual  variances,  and  correlation 
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coefficients  were  insensitive  to  a  geographic  grouping  of 
transformations.  The  difference  variances,  in  some  instances, 
rose  substantially  (Section  5.4).  This  was  confirmed  in  the 
spectra  at  Station  D  where  a  geographically  averaged  correction 
was  applied.  (Section  5.5) 

An  analysis  of  the  heat  flux  spectra  at  Stations  C  and  N 
demonstrated  seme  of  the  internal  properties  of  the 
transformations.  A  large  annual  cycle  in  both  the  latent  and 
sensible  3H  heat  fluxes  did  not  necessitate  a  seasonal 
transformation  at  Station  C  where  the  residual  variances  were 
lowest.  An  analysis  of  the  residual  spectra  at  Station  N, 
exactly  where  the  heat  flux  residuals  were  highest,  revealed 
that  a  seasonal  transformation  could  markedly  improve  the 
residual  variance.  Thus,  the  gross  fora  of  the  transformation 
may  be  dictated  by  the  vector  averaging  of  the  wind;  deviations 
from  the  broad  transformations  may  be  dependent  upon  the 
correlation  of  the  annual  air-sea  temperature  or  humidity 
differences  with  the  annual  winds.  Other  areas  of  high 
residuals  (eq.  Station  E,  latent  heat  flux)  may  also  require  a 
seasonally  dependent  transformation.  (Section  5.5) 

An  analysis  of  the  uncorrected  stress  component  spectra 
showed  a  consistent  underestimation  of  the  JH  spectra  by  the  VA 
spectra  at  all  frequencies.  Even  though  Station  D  showed  the 
highest  residuals  after  transformation,  there  was  no  evidence 
that  a  seasonal  transformation  improved  the  situation.  (Section 
5.5) 

Any  one  of  *he  three  +ransf or mations  ,T\  )  mill 
greatly  aid  calculations  based  on  air-sea  fluxes  calculated  from 
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averaged  constituent  data.  The  transf oraations  for  the 
constant  drag  coefficient  are  in  the  correct  range  to  aid  in 
improving  the  pressure  averaged  transports  of  Figure  1(b)  to 
■ore  closely  resemble  the  averaged  transports  of  Pigure  1(a). 
Consequently  application  of  the  transf oraation  aay  yield  tore 
accurate  results  from  numerical  models  wnile  using  lore 
realistic  air-sea  transfer  coefficients.  Furthermore,  the 


technique  offers  the  possibility  of  numerically  depicting  the 
time  history  development  of  specific  wind  forced  oceanic  events. 
Given  the  similarity  in  the  test  function  results  between 
,  and  \  •  the  two  empirical  formulae  should  ease 

implementation  of  the  transformations  with  little  loss  in 


accnracy . 

The  improvements  in  the  spectra  are  very  clear.  From  a 
situation  at  1=28.0  days  where  the  uncorrected  time  series 
yielded  large  discrepancies  between  the  VA  and  JH  spectra  and  a 
residual  spectrum  which  was  as  large  or  larger  than  the  VA 
estimated  spectrum,  the  )J  transformation  produced  a  match 
between  the  two  estimates  on  a  f requency-by-f reguency  basis  and 
produced  substantial  reductions  in  the  residual  power  density. 
(Section  5.5) 


Several  facets,  particularly  in  applying  transformations  to 
tar  erne trie  pressure  ca leu la t ions,  require  further  investigation. 


First,  the  barometric  pressure  maps  have  an  inherent  spatial 
averaging  which  could  not  be  simulated  in  this  study.  Similar 
calculations  at  the  eleven  weathership  locations  based  on 
qeostrophic  wind  data  (such  as  developed  by  Hi llebra n d ,  197 8) 
could  confirm  whether  similar  transformations  for  geostrophic 


151 


wind  are  required.  Application  of  the  transformations  to 
transport  calculations  similar  to  those  of  Aagaa rd  ( 1 970)  could 
add  further  confirmation.  Second,  the  sub-tropics,  tropics,  and 
mcnscon  regions  (where  a  strongly  deterministic  annual  signal  is 
inherent  in  the  winds)  have  not  been  investigated.  If 
confirmation  of  the  technique  could  be  gained  froa  replication 
on  geostrophic  winds,  then  it  may  be  possible  to  extrapolate 
transformations  to  other  regions  of  the  world  using  these  data. 

The  wind  rose  errors  were  investigated  at  all  levels  of 
construction.  In  all  cases,  the  average  errors  were  equal  to  or 
less  than  the  corrected  difference  mean  errors  for  the 
equivalent  L=28.0  days  of  the  vector  averaged  winds.  Slight 
overpredicticn  could  be  detected  in  the  discrete  versus  the 
continuous  calculation.  The  reduction  from  36  to  8  directions 
induced  a  sliqht  underestimation  in  the  3H  stress.  Errors  in 
the  MCA  (Marine  Climatic  Atlas)  calculation  were  probably 
affected  by  the  allocation  of  the  gale  force  winds  and  inherent 
correlation  between  the  winds  and  the  air-sea  temperature  and 
humidity  differences.  (Chapter  VI) 
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A^^endi x  A 


Development  of  some  mathematical  aspects  of  the  multiple 
regression  analysis. 


A .  1  Equivalence  Of  S  .An d  fl"' 


Let  X  he  the  3H  variate  and  X'  ne  the  VA  variate  and  k  be 
the  Beaufort  interval  (from  1  to  13).  Ihe  object  is  then  to 
minimize  the  residual  variance  defined  as: 

->  _  ‘2  r-  .  —  Z-  l 3  J .  1  ~7  2- 


h-(  i'l  **  *  UW--i M 


where  n  is  the  total  nuroter  of  points  in  the  averaging  length 
and  J  is  the  plumber  of  values  in  Beaufort  category  k. 

Setting  j.  =  Q  for  each  k: 


r-i  Ut  L  * 


Xa\'  ~  E  vte  X  n< 


)  =0  1. 


Redefine  : 

JV. 

^  h  -  X*  X*.^  X  *  ^ 
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r1  * 

T-  J- 

Xft  =  Zj  x«l\ 
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___  r*  , 

Xn  =  Zi 

1 3  Tk  v 
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Xi.  A* 


expanding  AX  in  full  gives: 


i3  r«. 


*  S|«.  -  Tk  S  Z!  Xi*^  -  £** 


1ms  equation  was  solved  by  Gaussian  elimination  with 
iterative  improvement. 

Expanding  equation  1.1A  gives: 


^  ^  (AxY 

^  la  ^  I  J 


Substituting  1.3A  gives: 


^  I  fv-  £  ♦  <a*y 


Similarly  the  difference  variance  is: 


(r«^=  i. 


Again  substituting  1.3A  gives: 

i-  •>  .  '  S  C  ■>  ^  t- - x. 


J  w £ lx* '  |j. +  «- zfn *) 

^  In.  ~  ^ 
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Or  : 


(Ty  AT, cl  - 


13 


2  *i.  t  Xft  (A*\ 

H  fe-.i  t 


it- 


-  (xixY 


1.  7  A 


'T  ATx  r  5T«  Vx  • 


Tnis  is  exactly  equal  to  1 .  5 A  or  ^ A^, '  <T»  ^ex  •  T*19 

approximation  for  the  wind  stress  was  that  all  linear  terms 
(means)  are  small  compared  to  the  quadratic  and  cross  terms  in 
which  case: 


I* 
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A.  2  Me  Effe c+  Of  The  Change  In  £e.e  Cn  Tx 

Let  X  and  X*  be  the  3H  and  WVA  variates,  andO+tV—  -I*  =  5.« 

*  v  ^  ™ 

,TxS»v«  is  tne  residual  variance  before  application  of 

T- 

andT  St>  is  the  residual  variance  after  application  of 

(Pfe')  .  Then : 


2.  1  A 


*  CTy  (T*(  -  3.  R.jt  G 1 T<  <Tv' 


2.  2A 


Where  is  the  averaging  operator. 

Which  reduces  to: 


;s 
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^  G  4  c')1"  Vi  ^ 


£  i + e)  <u- 
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Where  Exis  the  correlation  coefficient  =  (.xT1  -  H  k^/V*  1 
Since  Rv  is  invariant  under  a  scale  change  and 
(_  i*  fe')^.  ,  ^  ;  subtracting  2.3A  from  2.4A  gives 

Sxe')  -  +gt-_  2.  \^XG.C  »+<&)] 

=  £\>  +  €  -  ^6  +  e\j  r^« 

2.  5A 


Equation  2. 3A  can  be  rewritten  is  terms  of 


S  *  *  1.2*  2&  -  I’iZnCl+d)  £*"]  (f*( 


2.  6A 


Substituting  2.6A  into  2.  5A  gives 


<r. 


Dividing  by  <T*  £*<?  gives: 


C-S 


S\e 
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X  o 
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uustituting  2.5A  into  2.7A  gives: 

,  t -  f(e'*  v  e2-]’} 

\t  o 


2.8A 
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This  result  is  independent  of  whether  is  at  the 

minimum  or  not  and  thus  may  be  positive  or  negative. 

A.  3  The  Effect  Of  The  Change  In  E  On  when  'E  Kf  Is 

Minimized. 


Considering  the  primed  variates  to  be  uncorrected  and  tne 
unprimed  variates  to  be  the  varianles  corrected  for  the  DV  Dias, 
then  according  to  Section  A. 2: 

a-,1  *<r!  = 

3.  1  A 

^  ^  -  2  2k  <Tx 

At  the  minimum: 


f-e. 


■  < 


-  tr: 


*  T*  - 


which  reduces  to 


3.  2  A 


Now  : 

^  ^  y  -t -  rv  -  y. 


3.  3A 


<7^  (S  v<  -  Sve  }  -  v  "  V>  >  -2  R>(T2  (o*  -  o'*') 


Sucstit ut ing  Equation  3.1A  gives: 

--  -  C  -  *  (c rx-(r^  =  ({tv-r^ 


A.  4  j?rou_f  That  A  Shift  Of  l]/2  Is  Induced  Bet  wee  n  The  VA  And  JH 
Variates  At  Stati_cn  N». 


It  we  let  the  spectral  components  of  the  VA  series  at  any 
cycle  be  <a,*ifc,)  ,  and  the  3H  series  be  (a^+ibx).  Since  the 

courier  transform  of  the  residuals  is  a  linear  operation,  the 
power  at  the  residuals  is  : 

(a^a^  (j^ -\o^)  _  (a^ efcO'V i 6^2.4,  tO 

2_  2.  — - — - - - - —  4.  1 A 

7,  . 

if  we  consider  that  the  VA  variate  must  be  stretcned  and  rotated 
tc  match  the  2H  variate  then 

VCXv  +  H&,  +  ^  <9  4>2A 

where:  a.  m  is  the  stretch;  and 

b.  ^  is  the  angle  from  the  VA  tc  the  3H  spectral 

pairs. 
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Expanding  C  and  matching  real  and  imaginary  parts  gives: 
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Squaring  both  sides  of  Equation  4. 3A  gives  the  identity 


(at  -t-  bV)  =  yrf  (aN 


4.  4A 


Substituting  Equations  4. 2A  and  4. 3A  into  4.4A  gives 

(box-  0^Lj-  Ck-L-  '  Cc.,'1-1’  I' 


4.  5A 


From  Figure  5-12,  the  residual  power  at  the  annual  cycle  is 
approximately  equal  to  the  sum  of  the  pcwer  of  the  3H  and  VA 
series,  or: 

(o^-omT  .  (v-Vf* 

4.  6 A 


Using  4 . 4 A  gives: 

((XiVo^)  2  He  tos'©  =■  o 

Neither  m  nor  the  VA  power  density  equals  0.0  thus  9  =  tifjz 

or,  a  900  phase  shift  has  been  induced  between  the  3H  and  VA 
series  at  the  annual  cycle. 
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